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I.  Introduction 


This  report  only  discusses  work  that  has  been  carried  out  during  the  funding 
period  from  September  30,  1989  to  September  29,  1990  supported  under  Grant  No. 
AFOSR-85-0316.  This  grant  has  been  awarded  to  the  Polytechnic  University  since 
September  30,  1985.  The  work  accomplished  during  the  early  funding  period  has 
already  been  reported  in  the  previous  four  years'  progress  reports. 

The  investigation  covered  by  this  grant  deals  mainly  with  the  generation  of 
millimeter  waves  by  means  of  a  relativistic  electron  beam  and  with  the  interaction 
between  electromagnetic  waves  and  the  plasma.  However,  an  extension  of  the 
research  effort  to  study  the  mechanism  for  high  gradient  electron  beam  acceleration 
and  the  arc  discharge  in  high  current  carrying  plasma  has  also  been  made.  During 
the  funding  period  a  new  experimental  program  on  the  study  of  microwave 
absorption,  scattering  and  reflection  by  a  filamentary  plasma  produced  by  hollow 
cathode  discharge  and  frequency  shift  of  microwave  pulse  by  self -generated  plasma 
in  a  plexiglass  chamber  also  have  been  initiated.  Many  good  results  have  been 
obtained  and  will  be  presented  in  the  next  section. 

During  the  past  funding  year,  several  findings  of  our  work  have  already  been 
published  or  are  in  press.  Basically,  the  topics  of  investigation  are  divided  into  four 
general  categories:  (a)  optimization  of  the  design  of  a  compact,  high  power  mm 
wave  source  (cusptron  device)  and  its  operation,  (b)  generation  and  application  of 
artificial  atmospheric  plasma  including  propagation  of  high  power  microwave  pulse 
in  the  air  breakdown  environment  and  Bragg  scattering  study,  (c)  wave-plasma 
interaction,  study,  of  the  absorption,  scattering  and  reflection  of  microwaves  by 
hollow  cathode  discharge  plasma  and  (d)  frequency  up-conversion  of  microwave 
inside  a  suddenly  generated  plasma. 

For  the  first  topic,  particle  simulations  are  employed  in  the  analysis.  The 
study  focused  on  the  design  optimization  for  the  cusptron  device.  The  results  are 
then  used  for  the  experimental  setup  and  development  of  a  cusptron  (high 
harmonic)  microwave  tube.  This  tube  is  already  radiating  with  reasonable  output 
power  and  with  expected  output  frequency.  The  second  topic  is  devoted  to  the 
experimental  study  of  Bragg  scattering  of  EM  waves  by  microwave  produced  plasma 
layers  and  the  propagation  of  high  power  microwave  pulses  in  the  air  breakdown 


environment.  The  purpose  of  these  experiments  is  to  simulate  the  situation  of 
producing  plasma  layers  in  the  upper  atmosphere  for  over-the-horizon  (OTH) 
radars  application.  The  third  topic  includes  the  setup  of  a  hollow  cathode  discharge 
plasma  and  the  interaction  between  microwave  and  such  a  plasma.  The  last  one 
deals  with  an  initial  value  problem.  We  have  demonstrated  that  the  wave 
frequency  can  be  changed  in  a  time  varying  medium.  The  focus  of  study  is  on  the 
scattering  property  of  filamentary  plasma. 

Through  the  support  of  this  grant,  two  graduate  students  have  completed 
their  Ph.D.  degrees.  Dr.  K.  K.  Tiong  completed  the  requirements  of  the  degree  in 
August  1989.  His  Ph.D.  dissertation  was  included  in  the  last  year's  progress  report. 
Dr.  Y.  S.  Zhang  completed  the  requirements  in  August  1990.  His  Ph.D.  dissertation 
is  included  in  the  present  report. 
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Operational  characteristics  of  a  cusptron  microwave  oscillator 

K.  K.  TIONG,f  S.  P.  KUOt  and  P.  E.  MILLERf 

The  operational  cl  iracteristics  of  a  cusptron  oscillator  for  three  different  harmo¬ 
nic  interactions  (8,  10  and  12)  are  studied  via  numerical  particle  simulation.  The 
results  lead  to  the  determination  of  the  optimum  design  parameter  of  the  wave¬ 
guide  structure  for  both  TRNo,  and  TE02  modes.  For  the  TEWn)  modes,  it  is  found 
that  the  n  mode  is  roughly  twice  as  efficient  as  the  In  mode  counterparts.  The 
efficiency  study  of  the  tube  operating  at  the  TE02  mode  (where  only  the  2n  mode 
exists)  leads  to  a  single  optimum  waveguide  parameter,  b/a  =  1-3  (fo  and  a  are  the 
outer  and  inner  radii  of  the  waveguide)  for  all  three  harmonics. 


1.  Introduction 

Recently,  there  has  been  considerable  effort  in  generating  powerful  microwave 
radiations  by  relativistic  electron  beams.  Various  devices  including  gyrotrons  (Kuo 
et  al.  1986,  Granatstein  et  al.  1973),  relativistic  magnetrons  (Palevsky  and  Bekefi 
1979),  and  free-electron  lasers  (McMullin  and  Bekefi  1981)  have  been  studied  experi¬ 
mentally  and  theoretically.  In  general,  these  devices  are  not  aimed  to  be  physically 
compact  for  the  purposes  and  requirements  of  their  operation.  In  contrast,  cusp- 
trons  (Namkung  1984,  Kuo  et  al.  1988,  Tiong  et  al.  1987,  Lawson  and  Strifiler  1986. 
Chu  and  Dialetis  1984)  are  being  developed  with  the  aim  of  keeping  the  physical 
size  of  the  devices  compact. 

The  cusptron  tube  is  a  promising  device  for  efficient  operation  at  higher  cyclo¬ 
tron  harmonic  modes  (Namkung  1984,  Tiong  et  al.  1987,  Lawson  and  Strifiler  1986). 
A  cusp  magnetic  field  is  employed  to  produce  an  axis-encircling  electron  beam  (E- 
layer)  as  the  radiation  .mrce.  Such  a  beam  configuration  features  the  strongest 
finite  Larmour  radius  effect  (Chu  and  Dialetis  1984)  for  the  harmonic  interactions. 
This  device  also  uses  a  slotted  cylindrical  waveguide,  similar  to  the  anode  blocks  of 
a  magnetron,  to  determine  the  field  pattern  of  the  resonant  mode.  The  harmonic 
resonance  interaction  is  further  improved  by  the  periodic  structure  of  the  waveguide 
wall  which  introduces  periodic  fringe  fields  in  the  field  pattern  of  the  mode,  and 
enriches  the  resonant  harmonic  contents  of  the  mode  fields  as  experienced  by  the 
electrons  (Lau  and  Barnett  1982).  Moreover,  such  a  waveguide  structure  can  allow 
efficient  operation  even  at  low  beam  energy  and  minimize  the  competition  between 
the  modes. 

Similarly  to  the  magnetron  device,  the  cusptron  can  also  have  n  mode  and  2n 
mode  operations.  In  the  2 n  mode  operation,  the  phase  shift  from  one  slot  to  the 
next  is  360°,  and  thus,  the  number  of  slots  on  the  waveguide  wall  suggests  the 
optimum  cyclotron  harmonic  number  for  the  resonant  interaction.  In  comparison, 
the  phase  shift  between  two  adjacent  slots  is  180°  in  the  n  mode.  Hence,  n  mode 
operation  needs  twice  the  number  of  slots  for  the  same  cyclotron  harmonic  inter¬ 
action.  Consequently,  the  field  pattern  of  the  resonant  mode  in  each  operation  is 
expected  to  be  different  from  one  another.  This,  in  turn,  will  give  rise  to  a  difference 
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in  the  output  efficiency.  In  order  to  reach  a  conclusive  decision  on  the  mode  of 
operation,  a  comparative  study  of  the  performance  of  the  cusptron  tube  operating  in 
the  n  mode  and  2 n  mode  is  needed. 

In  the  present  work,  the  non-linear  evolution  of  the  radiation  fields  of  the  cusp¬ 
tron  tube  is  studied  by  using  a  particle  simulation  (Kuo  et  al.  1988).  The  functional 
dependence  of  the  efficiency  on  the  beam  energy  y,  and  on  the  ratio  b/a  of  the  outer 
radius  b  to  the  inner  radius  a  of  an  JV-slotted  waveguide  will  be  determined.  We  will 
consider  the  cases  of  N  =  8,  10  and  12  for  2n  mode  operation,  and  N  =  16,  20  and 
24  for  it  mode  operation.  The  purpose  of  this  study  is  to  identify  the  optimum 
design  parameters  and  the  mode  of  operation  for  the  device. 


X  Characteristic  equations 

The  resonant  interaction  between  an  axis-encircling  electron  beam  and  the 
modes  of  a  slotted  cylindrical  waveguide  is  studied.  The  cross-sectional  view  of  the 
waveguide  is  shown  in  Fig.  1.  The  structure  of  the  waveguide  is  characterized  by 
the  slot  number  N,  the  half-angle  of  the  slot  80 ,  and  the  radii  a  and  b.  Using  only 
the  fundamental  mode  fields  for  the  slot  regions,  the  boundary  conditions  of  the  TE 
mode  fields  at  p  =  a  lead  to  the  equation 


for  the  determination  of  the  transverse  wavenumber  k  of  the  mode  fields,  where  the 
cut-off  condition  kz  =  0  is  assumed;  T  =  m  4-  nN,  where  m  is  the  mode  index,  which 
determines  the  phase  difference  between  slots;  m  =  0  for  the  2 it  mode  and  m  =  N/2 
for  the  it  mode; 

.  =  [ Jo(kp)Y’0(kb)  -  Jo(kb)Yolkp)-] 

1  P)  \_J'0(ka)Y'0(kb)  -  J'0(kb)Y'0(ka)-\ 


Figure  1.  Cross-sectional  view  of  the  eight- vane  cylindrical  waveguide. 
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Figure  2.  Relationship  of  the  normalized  transverse  wavenumber  (ko)  with  the  waveguide 

parameter  ( b/a ). 


Je  and  y,  are  the  Bessel  and  Neumann  functions,  respectively;  a  50%  slot  opening  is 
considered  so  that  0O  =  n/2N. 

In  general,  (1)  has  multiple  roots  kj,j  =  1,  2,  3, ....  for  a  given  waveguide  struc¬ 
ture.  Each  root  is  a  function  of  the  ratio  b/a,  and  defines  a  frequency  range  of 
operation.  In  the  following  analysis  we  consider  for  each  N0th  (N0  =  N  for  the  2n 
mode  and  N0  =  N/2  for  the  n  mode)  harmonic  interaction,  three  branches  of  roots 
starting  from  b/a  =  1,  namely  the  TE*ol  (both  7t  and  2jt  modes)  and  TE02  (2n  mode 
only)  modes  respectively.  These  branches  of  roots  are  presented  in  Fig.  2. 

To  achieve  the  resonant  interaction,  we  have  to  adjust  the  background  magnetic 
field  to  the  strength  such  that  ft  &  to0/N0  for  the  N0th  harmonic  cyclotron  reson¬ 
ance,  where  ft  =  fto/y  is  the  relativistic  electron  cyclotron  frequency,  ft0  =  eB0/mc, 
and  o)0  is  the  eigenfrequency  of  the  empty  waveguide  mode.  The  motion  of  electrons 
in  the  wavefields  and  a  uniform  axial  magnetic  field  B0t  is  governed  by  a  set  of 
coupled  non-linear  equations 


Z  R  =  P/ym 
at 

(2) 

JP=I 

dr 

|e  +  J»x(B  +  f  B0)J 

(3) 

where  y  =  (1  +  P2/m2c2)112,  P  =  ymv,  and  E  and  B  are  the  waveguide  TE  mode 
fields. 

This  set  of  equations  can  be  simplified  only  if  the  resonant  component  of  the 
total  force  is  included  on  the  right-hand  side  of  (3).  This  approximation  is  justifiable 
since  the  perturbations  on  the  electron  orbits  are  mainly  introduced  by  the  resonant 
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components  of  the  total  force  which  have  the  slowest  time  variation.  Consequently, 
only  the  resonant  trajectory  ( R ,  <f>)  will  be  of  interest.  This  trajectory  can  be 
expressed  to  be 


r  =  RL(t)  =  w«n0 
<!>  =  <t>  o  +  I m  +  fW)  df' 


where  RL(t)  and  \p(t),  to  be  determined  self-consistently,  are  the  time-dependent 
Larmour  radius  and  the  slowly  varying  phase  of  the  electrons:  this  results  from  the 
resonant  interaction  between  the  electrons  and  wave  fields.  From  (2)  and  (3),  the  set 
of  self-consistent  characteristic  equations  (Chu  and  Dialetis  1984,  McMullin  and 
Bekefi  1981)  for  the  resonant  trajectory  of  the  electrons  is  derived  to  be 


^  a  =  AE0J'No(a)  cos  <t>  (4) 

d  f  Aco  Id  1 

—  <£  =  -  Aco  +  AE0< —  J N  (a)  =  -—  [aJ'No(a)]  \  sin  <t>  (5) 

at  [  a)  a  da  ) 


where  a  =  kKL(t)  and 


4>  =  N0(*  +  *0)- 


i 


AaAjt')  dt' 


£0  is  the  amplitude  of  the  wavefield; 

Act)  =  a)  —  N0tl  =  Act)0  +  o>o(y  ~  7o/y) 


Act)0  a*  cu  —  a)0 ,  and 

o>o  ~  N0Q0/y0;  A  —  aNow/B0 ,  aNo  =  1  /N0JNo(ka) 

and  y  and  a  are  related  through  the  relation 

y2  -  l  +  (fWco)2 


The  z  component  of  (3)  reduces  to  (d/d t)Vz  =  0  for  kz  =  0  and  is  not  coupled  to  the 
rest  of  the  components.  A  similar  result  is  also  obtained  under  the  grazing  condition 
that  vz  =  k,  c2/a>  =  vt .  The  effect  of  finite  kt  and  vz0  on  optimum  interaction  is 
found  to  be  small  as  long  as  kxc  <  to. 

Since  the  field  amplitude  £0  varies  during  the  interaction  with  electrons,  it  has  to 
be  determined  self-consistently.  This  consideration  can  be  included  by  coupling  (4) 
and  (5)  to  the  energy  conservation  equation,  which  is  given  by 


(6) 


where  <  >  represents  an  average  over  the  initial  random  phase  distribution  of  elec¬ 
trons;  n0  is  the  electron  beam  density  averaged  over  the  inner  cross-section  of  the 
waveguide  structure  and  defines  an  averaged  plasma  frequency 


(Up  =  (4 nn0e2/m)t/2 


Cusptron  microwave  oscillator 


131 


€  is  the  effective  dielectric  constant  of  the  waveguide  structure  for  the  TE  mode  and 
is  given  by 

sin  rfl0Tr  _1_  Mkaj  (  Jfjkaf  | 

T  J  |_  ka  J'xika)  \  k2a2)  J'^(ka)\ 

+  £2  [£  G‘(kh)  -  GH ka)  -  (1  -  SJ  (7) 

with 


G(kp)  =  F"(kp)  - 


J\lkp)Y,(kb)  -  J,(kb)Y\(kp) 
J\(ka)Y1(kb)  -  J^kb^ka)  ’ 


a  =  b 
a  ¥=b 


Equations  (4)  and  (5)  describe  the  motion  of  a  single  electron  in  the  wave  fields, 
and  (6)  describes  the  evolution  of  the  self-consistent  field  driven  by  the  collective 
response  of  the  electrons.  In  order  to  couple  them  together,  a  phase  average  over  the 
single  electron  motion  is  necessary.  This  will  be  done  numerically,  in  the  following 
section,  by  solving  the  trajectories  of  a  large  number  of  electrons  with  random  initial 
phase  distribution. 


3.  Numerical  simulation 

Initially,  the  particles  are  set  to  be  positioned  symmetrically  about  the  axis  of  the 
waveguide  structure,  and  the  energy  of  each  particle  is  given  by  a  predetermined 
value  resulting  from  the  cathode  voltage.  This  distribution  represents  an  unper¬ 
turbed  sample  of  the  actual  randomly  distributed  electrons  in  the  ideal  £-layer 
which  have  an  initial  density  per  unit  length  of  NL  =  1-6  x  IO10  cm-1. 

With  each  increment  of  simulation  time  during  program  execution,  the  pertur¬ 
bation  in  each  particle’s  position  is  computed  by  integrating  the  coupled  equations 
(4)  and  (5);  from  this  calculation,  the  particle  energy  can  also  be  determined.  By 
comparing  the  resulting  average  particle  energy  with  the  initial  average  particle 
energy,  the  energy  transferred  to  the  wavefield  can  be  determined  through  (6).  The 
next  time  increment  repeats  the  calculations  using  the  perturbed  position  and 
energy  values  for  the  particles,  along  with  the  resulting  wavefield  magnitude. 

In  this  simulation,  it  was  found  that  128  particles  were  sufficient  to  obtain  accu¬ 
rate  data;  greater  numbers  of  particles  only  resulted  in  longer  execution  time.  The 
program  was  run  on  a  CRAY  X-MP/48  operating  at  the  Pittsburgh  Supercomputer 
Center.  The  coupled  equations  (4)  and  (5)  were  solved  using  the  ODEPACK 
LSODE  Fortran  subroutine  available  on  the  CRAY. 


4.  Numerical  results 

In  this  section,  we  study  the  operational  characteristics  of  cusptron  devices  by 
using  a  computer  particle  simulation  approach  described  in  the  previous  section.  In 
all  the  calculations,  we  are  using  an  ideal  £-layer  having  a  constant  line  density  of 
the  typical  value  NL  =  7ta2n0  =  1-6  x  lO^cm'1  as  the  radiation  source.  However, 
the  initial  beam  energy  (y0)  will  be  considered  as  a  variable.  In  addition,  the  ratio  of 
the  outer  to  inner  radius  of  the  waveguide,  b/a,  is  also  varied  in  the  numerical 
analyses  for  three  harmonic  cases  ( N0  —  8,  10  and  12).  As  shown  in  Fig.  2,  the  mode 
frequency  w0  ~  kc  varies  with  b/a,  so  does  the  cyclotron  frequency  fl0  - 
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(a) 


(b) 

Figure  3.  (a)  Efficiency  versus  b/a  for  the  In  mode  (N0  =  8(0)  operating  at  the  TEWo,  mode 
with  normalized  Larmour  radius  RJa  =  0-7,  originated  from  J'N  =  0,  first  zero.  (b). 
Efficiency  versus  b/a  for  the  n  mode  {N0  =  8,  10  and  12;  N  =  16  (□),  20  (A)  and  24 
(O))  operating  at  the  TEWo,  mode  with  normalized  Larmour  radius  RJa  =  0-7,  orig¬ 
inated  from  J'h, 2  =  0,  first  zero. 
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y0  (o0/N0  ~  y0  kc/Na .  In  the  experimental  arrangement,  the  radius  of  the  £-layer 
produced  by  the  cusp  magnetic  field  is,  in  general,  fixed  by  the  radius  of  the  ring 
cathode  of  a  pierce  gun.  Therefore,  the  beam  energy  (y0)  has  to  vary  with  b/a  in 
order  to  match  the  electron  Larmour  radius  RL  to  the  radius  R  of  the  cathode, 
otherwise,  the  produced  axis-encircling  beam  may  either  become  non-concentric 
with  large  spatial  spread  (for  y0  to  be  too  low)  or  have  large  vz0  (for  y0  to  be  too 
large).  In  either  case,  the  performance  of  the  tube  will  be  degraded.  In  order  to 
optimize  the  interaction,  one  would  like  to  set  the  radius  R  of  the  cathode  as  close 
as  possible  to  the  value  a  of  the  inner  radius  of  the  waveguide.  This  is  because  the 
high-harmonic  mode  fields  are  more  concentrated  in  the  region  near  the  inner  wall 
of  the  waveguide.  However,  for  practical  application,  one  would  have  to  keep  a  gap 
to  allow  for  the  possible  spread  of  the  E-layer.  The  typical  value  of  RJa  is  chosen  to 
be  between  0-7  and  0-8.  For  this  reason,  the  beam  energy  used  in  the  calculation  will 
be  varied  with  b/a  to  ensure  that  0-7  <  RJa  ^  0-8. 

Presented  in  Figs  3(a)  and  ( b )  are  the  functional  dependences  of  the  saturation 
efficiency  rj  on  the  ratio  b/a  for  the  Jrth  mode  and  2 n  mode,  respectively.  In  each 
mode  of  operation,  three  harmonic  cases  (N0  =  8,  10  and  12)  are  considered,  where 
the  beam  radius  is  fixed  to  be  0-7 a  and  k  is  chosen  from  the  upper  branches  of  the 
dispersion  curves  of  Fig.  2,  whose  k  values  at  b/a  =  1  correspond  to  those  of  TENol 
modes.  The  results  in  Fig.  3(a)  can  be  used  to  determine  the  most  favourable  b/a 
value(s)  for  each  N0  from  which  the  optimum  output  efficiency  can  be  achieved  for 
2n  mode  operation.  Those  values  of  b/a  are  identified  to  be  b/a  =  2  0  for  the  N0  =  8 
case,  b/a  =  1-5  for  N0  =  10,  and  b/a  =  1-3  and  1-9  for  N0  =  12. 

Using  these  b/a  values  for  the  2n  mode  tubes,  their  operating  frequencies  are 
defined.  In  order  to  make  a  meaningful  comparison  between  the  performance  of  the 
n  mode  and  2 n  mode,  their  operating  frequencies  should  be  set  to  be  the  same  for 
each  harmonic  case.  The  inner  radius  a  of  the  waveguides  and,  hence,  the  beam 
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energy  are  also  kept  the  same  for  both  modes  of  operation.  Thus,  the  corresponding 
b/a  values  for  n  mode  operation  can  be  determined  from  the  dispersion  curves  of  the 
it  mode  with  the  ks  of  the  2 n  mode.  Figure  4  displays  the  functional  dependence  of 
the  saturation  efficiencies  q  on  the  beam  energy  y0  for  the  specified  b/a  values,  where 
y0  only  varies  in  the  range  for  which  0-7  ^  RJa  <  0-8.  It  shows  that  the  n  mode 
always  has  a  better  output  efficiency  than  that  of  the  2 n  mode. 

The  frequencies  defined  by  the  lower  branches  of  the  curves  of  Fig.  2  correspond 
to  the  modes  branching  out  of  the  TE02  mode  due  to  the  presence  of  slots  on  the 
wall.  Thus,  tubes  operate  only  at  27t  mode  for  these  dispersion  branches.  Again,  we 
first  determine  the  dependence  of  r\  on  b/a  for  each  N0  case.  The  results  are  present¬ 
ed  in  Fig.  5  for  RJa  =  0-7.  The  most  favourable  b/a  values  for  the  design  of  tubes 
can  then  be  obtained  from  the  locations  of  the  peaks  of  those  curves.  These  values 
are  found  to  be  quite  insensitive  to  N0  and  are  all  about  equal  to  1-3.  This  is  shown 
by  the  fact  that  all  the  dispersion  curves  are  very  close  to  each  other.  Thus,  the 
harmonic  content  introduced  by  the  slots  depends  mainly  on  b/a ,  but  not  the 
number  of  the  slots.  Using  this  common  value  for  all  three  N0  cases,  the  dependence 
of  the  optimum  efficiency  on  y0 ,  which  is  restricted  to  vary  in  the  range  for  which 
0-7  ^  RJa  ^  0-8,  is  determined  for  each  N0 .  The  results  are  presented  in  Fig.  6. 

5.  Summary 

We  have  used  a  numerical  particle  simulation  to  study  the  non-linear  evolution 
of  the  radiation  fields  excited  through  the  negative  mass  instability  of  an  ideal  £- 
layer,  which  is  rotating  inside  a  slotted  waveguide.  The  dependences  of  the  satura¬ 
tion  efficiency  of  the  cusptron  device  on  the  structure  of  the  waveguide  (i.e.  b/a  and 
N0),  the  mode  of  operation  (n  mode  or  2 n  mode),  and  the  electron  beam  energy  (y0) 


Figure  5.  Efficiency  versus  b/a  for  the  2n  mode  (N0  =  8,  10  and  12)  operating  at  the  TE02 
mode  with  normalized  Larmour  radius  RJa  -  01.  Originated  from  J'0  =  0,  second 
zero. 


Cusptron  microwave  oscillator 


135 


Figure  6.  Efficiency  versus  relativistic  factor  (y)  for  the  2n  mode  operating  at  the  TE02 
mode.  0-7  <  RJa  <  0-8,  N  =  8  (□),  10  (A)  and  12  (<>)■  Originated  from  J'0  =  0, 
second  zero. 

are  thus  determined.  Three  branches  of  waveguide  modes  whose  dispersion  relations 
are  shown  in  Fig.  2  have  been  considered.  A  comparison  between  the  performance 
of  n  mode  and  2 n  mode  operation  has  also  been  made. 

The  present  study  leads  to  the  following  conclusions. 

(a)  Operating  the  waveguide  mode  in  the  harmonic  branch  (i.e.  it  reduces  to  the 
TENo1  mode  for  b/a  =  1),  the  jt  mode  tube  (N  =  2N0)  is  found  to  be  about  twice 
as  efficient  as  the  2ji  mode  tube  ( N  =  N0). 

(b)  The  optimum  structures  of  the  2ji  mode  tubes  are  determined  to  be  b/a  =  2  for 
N0  =  8,  b/a  =  1-5  for  N0  =  10  and  b/a  =  1-3  for  N0  =  12.  The  optimum  effi¬ 
ciency  for  each  corresponding  tube  is  found  to  be  15-5%  with  y0  =  114,  5-5% 
with  y0  =  1-36,  and  3%  with  y0  =  1-75. 

(c)  In  the  fundamental  branch  (i.e.  the  waveguide  mode  becomes  TE02  for  b/a  =  1), 
tubes  can  only  operate  at  2 it  mode.  A  common  value  of  b/a  =  1-3  is  found  to  be 
the  optimized  tube  structure  for  all  three  N0  cases.  The  optimum  efficiency  is 
found  to  be  about  18%  for  all  three  different  N0  tubes,  except  the  required  beam 
energy  for  each  tube  is  different,  y0  =  1-15  for  N0  =  8,  y0  =  1-09  for  N0  =  10 
and  y0  =  1-065  for  N0  =  12. 
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Operation  of  a  high  harmonic  cusptron  oscillator 

K.  K.  TIONGf  and  S.  P.  KUOf 

A  cusptron  oscillator  has  been  set  up.  It  uses  a  16-slot  waveguide  for  the 
generation  of  the  16th  cyclotron  harmonic  TE02  mode  (2n-mode  operation)  and 
the  8th  cyclotron  harmonic  TESI  mode  (jt-mode  operation).  The  output 
frequencies  are  12-5  GHz  and  8-33  GHz,  respectively.  Moreover,  a  15th  cyclotron 
harmonic  TE0J  mode  has  also  been  generated  at  10 GHz  via  peniotron 
instability. 

1.  Introduction 

It  is  evident  from  most  experimental  indications  that  gyrotrons  generally 
produce  good  operating  efficiencies  in  the  fundamental  and  second  harmonic 
modes  (Carmel  et  al.  1983,  Kreischer  and  Temkin  1987,  Granatstein  and  AlexefT 
1987).  However,  when  it  is  operating  in  higher  harmonic  modes  or  the  whispering 
galley  modes,  the  efficiency  drops  dramatically.  The  information  on  the  saturation 
efficiency  of  gyrotron  operation  at  higher-order  modes  is  therefore  of  considerable 
interest  to  a  number  of  researchers  (Chu  1978,  McDermott  et  al.  1983,  Kuo  et  al. 
1986).  Among  them,  Kuo  et  al.  (1986)  performed  a  non-linear  analyses  of  gyrotron 
operation  at  harmonics  as  high  as  N  =  8.  The  results  showed  a  progressive 
decrease  in  efficiency  with  increasing  harmonic  number.  However  the  efficiency 
increases  with  increasing  beam  energy.  The  demand  for  high-power  high-frequency 
sources  motivates  the  search  for  many  other  novel  ways  of  improving  the 
operation  of  high  harmonic  gyrotrons.  Further,  it  is  also  desirable  that  new 
devices  should  be  reasonably  compact,  for  applications  such  as  space  based  radar 
systems,  and  also  should  offer  a  degree  of  tunability  either  in  the  CW  range  or  the 
step  range.  Destler  et  al.  (1981)  showed  that  the  efficiency  of  a  gyrotron  operating 
at  higher  harmonics  can  be  improved  by  using  a  high  energy  rotating  beam.  The 
improvement  lies  in  the  fact  that  the  Larmor  radius  of  the  electron  beam  is 
directly  proportional  to  its  energy.  Therefore,  by  increasing  the  energy,  the  electron 
beam  will  be  brought  closer  to  the  wall  where  the  harmonic  field  is  concentrated 
(Lawson  et  al.  1985).  However,  this  feature  will  make  the  device  rather  bulky  when 
the  accelerator  is  required  to  produce  electrons  of  energy  higher  than  500  keV. 
Other  complications  also  arise  from  using  high-voltage  high-current  beams.  For 
example,  the  cavity  size  must  be  scaled  correspondingly  to  handle  the  increase  in 
dissipated  power  and  high  current  operation  in  overmoded  cavities  leads  to  the 
problem  of  the  suppression  of  parasitic  modes  which  would  otherwise  complicate 
the  mode  purity  controllability. 

To  improve  the  efficiency  of  a  high  harmonic  gyrotron  without  using  a  high 
energy  beam,  a  cusp-injected  hollow  beam  system  is  introduced.  The  system  has 
the  attraction  of  maximizing  the  finite  Larmor  radius  effect,  which  will  enhance  the 
beam-high-order  modes  interaction  without  using  a  high  energy  beam.  Such  a 
beam  trajectory  can  be  produced  by  passing  an  axuaily  streaming,  thin  annular 
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electron  ring  through  a  cusp-magnetic  field  configuration  created  by  magnetic  field 
reversal.  Using  this  system  and  with  an  appropriate  guiding  magnetic  field  in  the 
interaction  region,  the  E-layer  produced  can  be  positioned  as  close  to  the 
waveguide  wall  as  physically  permissible.  However,  other  disadvantages  including 
the  loss  of  electrons  in  traversing  such  a  field  configuration  and  radial  beam 
broadening  may  also  negate  some  of  the  benefits  of  this  system.  It  can  also  be 
shown  that  the  selectivity  of  the  higher-mode  operation  can  be  improved  by  the 
inclusion  of  periodic  boundary  slots  in  the  resonant  circuit.  As  the  ratio  of  b/a 
exceeds  a  certain  value,  which  depends  on  the  number  of  slots,  the  slotted 
boundary  creates  a  large  periodic  fringe  field  near  the  orbits  of  the  concentrically 
circulating  electrons,  thereby  enriching  the  harmonic  contents  of  the  RF  fields  as 
experienced  by  the  electrons  (Lau  and  Barnett  1982,  Chu  and  Dialetis  1984, 
Lawson  and  Strifller  1986,  Tiong  et  al.  1988).  These  added  features  have  the 
advantage  of  making  the  device  physically  compact  as  the  magnetic  field  strength 
required  for  high-harmonic  radiation  will  be  decreased  proportionally  to  the 
harmonic  number  and  the  required  beam  energy  for  a  given  Larmor  radius  will 
also  be  reduced  by  the  square  of  the  harmonic  number. 

The  experimental  feasibility  of  such  a  device  incorporating  the  previously 
described  features  was  first  demonstrated  by  NamKung  (1984).  In  his  work,  a  six- 
vane  resonant  circuit  with  b/a  =1-5  was  used.  Radiation  at  4-38  GHz  with  500  W 
power  was  detected  with  a  magnetic  field  strength  of  260  G.  In  his  later  work,  he  is 
able  to  generate  radiation  at  10-4  kW  power  at  6  GHz  (NamKung  et  al.  1987).  The 
magnetic  field  strength  required  was  350  G.  In  the  later  report,  it  was  noted  that 
the  4-38GHz  radiation  corresponds  to  the  4th  harmonic  and  the  efficiency  is 
optimized  to  be  9%  at  the  9kW  power  level,  while  that  of  6  GHz  corresponds  to 
the  6th  harmonic  radiation.  Later  reports  by  Chojnacki  et  al.  (1987),  using  the 
same  system  indicated  radiation  at  6-78  GHz  at  the  8th  harmonic  and  8-85  GHz  at 
the  10th  harmonic,  where  the  radiated  powers  were  reported  as  9-1  kW  and 
2-2  k W,  respectively. 

In  the  present  work,  we  extend  the  cusptron  operation  to  even  higher 
harmonics.  It  is  noted  that  in  all  previous  work  on  the  cusptron,  the  mode  of 
operation  is  on  the  TE0,  branch  of  the  dispersion  curve.  It  is  also  our  aim  to 
operate  the  device  in  the  TEN1  or  TE0M  branches  of  the  dispersion  curves  (i.e.  the 
solution  of  the  waveguide  mode  originating  from  the  first  zero  of  the  Nth 
harmonic  or  Mth  zero  of  the  fundamental  mode).  The  main  advantage  of  such 
operation  is  the  resultant  high-frequency  output  with  lower  magnetic  requirement. 
Other  objectives  include  exploring  the  operation  of  the  device  in  the  7i-mode 
configuration,  which  has  been  shown  by  particle  simulation  (Kuo  et  al.  1988, 
Tiong  et  al.  1990)  to  be  more  efficient  in  comparison  to  the  more  common  271- 
mode  operation. 

In  §2,  descriptions  of  the  background  theories  and  particle  simulation  for  the 
design  of  the  tube  are  presented.  The  experimental  set-up  and  results  are  presented 
in  §3.  Final  conclusions  and  some  perspectives  of  the  work  are  discussed  in  §4. 


2.  Characteristic  equations  and  particle  simulations 

The  transverse  wave  number  k  of  the  TE  mode  fields  in  an  N-slot  cylindrical 
waveguide  is  given  by  (Kuo  et  al.  1988) 
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where  the  cut-off  condition  fc,=0  is  assumed;  T  =  m  +  niV;  m  is  the  mode  index 
which  determines  the  phase  difference  between  slots;  m  =  0  for  2 it  mode  and  m  =  N/ 
2  for  it  mode;  90  is  the  half  angle  subtended  by  each  slot; 

F(k) = :j0(kp)r0(kb)-Jo(kb)Y0(kpn 
P)  lf0(ka)Y'0(kb)-J'0(kb)Y'0(kafi 

a  and  b  are  the  inner  and  outer  radii  of  the  waveguide;  J,  and  Y,  are  the  Bessel 
and  Neumann  functions,  respectively;  a  50%  slot  opening  is  considered  so  that 
90  =  it/2N. 

In  general,  (1)  has  multiple  roots  k/,  j=  1,2,3,...  for  a  given  waveguide 
structure.  Each  root  is  a  function  of  the  ratio  b/a,  and  defines  a  frequency  range  of 
operation.  The  dashed  curves  I  and  II  presented  in  Fig.  1  represent  the  roots 
corresponding  to  the  TE02  dispersion  branch  (2n  mode  only)  and  the  TE81 
dispersion  branch  {it  mode  only),  respectively.  For  the  it  mode  case,  only  the  range 
1  <b/a<  1-4  is  of  interest  since  ka  decreases  drastically  with  b/a. 

To  achieve  resonant  interaction,  we  have  to  adjust  the  background  magnetic 
field  to  a  strength  such  that  S}=a>0/N0  for  the  N0th  harmonic  cyclotron 
resonance,  where  Q=Q0/y  is  the  relativistic  electron  cyclotron  frequency  and  a>0  is 
the  eigenfrequency  of  the  empty  waveguide  mode.  As  the  resonance  condition  is 
met,  the  perturbations  on  the  electron  orbits  are  mainly  introduced  by  the 
resonant  component  of  the  total  force  which  has  the  slowest  time  variation. 
Consequently,  only  the  resonant  trajectory  will  be  of  interest.  A  set  of  self- 
consistent  characteristic  equations  for  the  resonant  trajectory  of  the  electrons  is 
thus  derived  to  be  (Kuo  et  al.  1988) 

a  =  AE0  J'jy0(a)  cos  <p  (2) 

y  <p= -Acu-ME0{(Ao)/a))J* 0(a)+-  [aJ'Wo(a)]}  sin  <p  (3) 

at  a  dot 


where 


a  —  kvJQ.  and 


<P  =  No('l'  +  <t>o)~ 


i: 


Ac o(t)  dt 


E0  is  the  amplitude  of  the  wavefield; 


Aco  =  co  —  N0Q  =  Aa>0  +  a»0(y  -  y0)/v,  Acu0  =  co  -  co0,  and  co0  =  NCi0/y0 

and  y  and  a  are  related  through  the  relation  y2=  1  +(Ci0ot/co)2. 

Since  the  field  amplitude  E0  is  going  to  vary  during  the  interaction  with 
electrons,  it  has  to  be  determined  self-consistently.  This  consideration  can  be 
included  by  coupling  (2)  and  (3)  to  the  energy  conservation  equation,  which  is 
given  by 
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(e/8 na2N0)  E%  =  n0mc2  <y> 


dt 


dt 


(4) 


where  *< )’  indicates  an  average  over  the  initial  random  phase  distribution  of 
elections;  n0  is  the  electron  beam  density  averaged  over  the  inner  cross-section  of 
the  waveguide  structure;  e  is  the  effective  dielectric  constant  of  the  waveguide 
structure  for  the  mode  under  consideration  and  is  given  by  (Kuo  et  al.  1988) 


.=  y  /sinrfloy  f  \Jr(ka)  (  T2  \J2(ka)} 
'  tJzm\  r  j  \  kaJ’r(ka)  V  k2a2J  Jr(ka)j 


+ 


2 N2 


G\kb)-G2(ka ) 


(5) 


with  G(kp)  —  F“(kp)  and  bah  is  the  Kronecker  delta. 

In  order  to  couple  (2)— (4)  together,  a  phase  average  over  the  single  electron 
motion  is  necessary.  This  is  achieved  numerically  by  solving  the  trajectories  of  a 
large  number  of  electrons  with  a  random  initial  phase  distribution.  Initially,  the 
particles  are  set  to  be  positioned  symmetrically  about  the  axis  of  the  waveguide 
structure,  and  the  energy  of  each  particle  is  given  by  a  predetermined  value 
resulting  from  the  cathode  voltage.  This  distribution  represents  an  unperturbed 
sample  of  the  actual  randomly  distributed  electrons  in  the  ideal  E-layer.  With  each 
increment  of  simulation  time  during  program  execution,  the  perturbation  in  each 
particle’s  position  is  computed  by  integrating  (2)  and  (3);  from  this  calculation,  the 
particle  energy  can  also  be  determined.  By  comparing  the  resulting  average 
particle  energy  to  the  previous  average  particle  energy,  the  energy  transferred  to 
the  wavefleld  can  be  determined  through  (4).  The  next  time-increment  repeats  the 
calculations  using  the  perturbed  position  and  energy  values  for  the  particles,  along 
with  the  updated  wavefield  amplitude.  In  this  simulation,  it  was  found  that  128 
particles  were  sufficient  to  obtain  accurate  data.  Adding  more  particles  in  the 
simulation  will  not  improve  the  results  over  the  accuracy  of  the  numerical 
scheme — Los  Alamos  ordinary  differential  equation  solver  (LSODE).  The  temporal 
evolution  of  the  wavefield  amplitude  can  thus  be  determined  and  the  saturation 
efficiency  t]  of  the  output  radiation  can  be  calculated.  The  solid  curves  I  and  II  in 
Fig.  1  show  the  functional  dependence  of  rj  on  b/a  for  27t-mode  16th  harmonic 
operation  (N0=16)  in  the  TE02  mode,  and  for  r-mode  8th  harmonic  operation 
(N0  =  8)  in  the  TEgl  mode,  respectively.  In  the  calculation,  the  beam  energy 
changes  with  b/a  so  that  the  ratio  rja  can  be  fixed  at  a  constant  value  0-7,  where 
rL  is  the  electron  Larmor  radius.  This  shows  that  the  first  peak  of  the  efficiency 
(solid)  curve  I  is  located  near  b/a  —  1-3.  The  efficiency  curve  II  for  the  TEgI  n  mode 
is  a  monotonicaliy  increasing  function  of  b/a  with  severe  trade-off  on  the  output 
frequency  at  b/a  =  1-3;  the  efficiency  of  the  n  mode  is  already  much  higher  than 
that  of  the  2 n  mode.  On  the  other  hand,  the  operation  frequency  of  n  mode  is  still 
comparable  to  that  of  2jt  mode.  Therefore,  b/a=  1-3  is  considered  to  be  a  good 
parameter  for  the  design  of  a  slotted  waveguide  able  to  be  used  for  In-  and  n- 
mode  operation  with  the  same  experimental  set-up.  The  second  peak  of  efficiency 
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b/a 

Figure  1.  Efficiency  and  ka  versus  b/a  for  the  2n  and  rr  inodes. 


curve  I  is  not  considered  because  it  requires  a  larger  b/a  value  and  also  operates  at 
a  much  lower  output  frequency  for  both  the  n  and  In  modes. 


3.  Experiments  and  results 

In  the  experimental  set-up,  twelve  individually  wound  copper  wire  coils 
connected  to  three  power  supplies  are  used  to  produce  the  desired  magnetic  field 
profile.  A  cusp  field  centred  at  the  anode  location  is  produced  by  the  three  coils 
positioned  in  the  upstream  region  together  with  the  first  three  of  the  remaining 
nine  placed  in  the  downstream  region.  The  other  six  coils  in  the  downstream 
region  are  used  to  produce  a  uniform  field  inside  the  waveguide  for  cyclotron 
harmonic  resonance  interaction.  Each  coil  can  sustain  a  maximum  of  20  A  current 
(~01T  magnetic  field)  without  overheating  for  an  extended  period  of  operation. 
The  cusp  transition  width  is  substantially  narrowed  by  a  soft  iron  plate  positioned 
between  the  upstream  and  downstream  regions  and  jointed  as  part  of  the  anode. 
The  cusp  transition  width  is  measured  to  be  about  0-5 cm.  The  system  vacuum  is 
maintained  by  an  ion  pump  at  about  IO~8Torr.  A  schematic  of  the  set-up  is 
shown  in  Fig.  2. 

As  discussed  in  §2,  b/a  =  1-3  is  a  good  choice  for  a  16-slot  waveguide.  The 
value  b/a  =1-3  is  therefore  picked  for  the  construction  of  the  10cm  long  slotted 
waveguide  whose  photo  is  shown  as  Fig.  3.  The  corresponding  ka  values  found 
from  the  dashed  curves  I  and  II  of  Fig.  1  are  5  14  for  the  TE02  In  mode  and  3-43 
for  the  TEgl  7r  mode,  respectively.  The  inner  radius  a=  l -98  cm  is  chosen  so  that 
0-65 <rja< 0-75,  where  rc  is  the  radius  of  the  annular  cathode.  A  hollow  electron 
beam  is  produced  from  an  annular  thermionic  cathode  of  radius  1-4  cm  and  radial 
thickness  0-16cm  with  a  Pierce-type  focusing  electrode.  The  electron  axial  energy 
is  provided  by  a  hard  tube  modulator  with  a  maximum  capability  of  30  kV  and 
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Figure  2.  Schematic  of  the  cusptron  experiment. 


Figure  3.  Photograph  of  16-slot  waveguide  (b/a=  1-3). 

10  A  current.  The  pulse  width  and  pulse  repetition  rate  are  adjustable  from  10-2 /rs 
and  200  Hz-1  kHz,  respectively.  The  cathode  is  mounted  on  bellows  coupled 
tubing  to  facilitate  alignment  and  the  adjustment  of  the  cathode-anode  gap 
without  breaking  the  vacuum  system.  The  anode  with  an  annular  slit  of  width 
0-2  cm  supported  by  three  bridges  is  attached  to  the  iron  plate.  The  slit  will  allow 
the  hollow  beam  to  pass  through  the  magnetic  cusp  transition  region,  where  the 
v.x  B  /orce  converts  the  beam  axial  velocity  v ,  into  the  azimuthal  velocity  t>,  on 


High  harmonic  cusptron  oscillator 


7 


the  downstream  side  of  the  cusp  transition.  The  16- vane  waveguide  for  the  beam- 
wave  interaction  is  placed  immediately  behind  the  anode  in  the  downstream 
region. 

Since  the  radius  rL  of  an  ideal  E-layer  is  equal  to  the  radius  of  the  hollow 
beam,  it  is  thus  fixed  by  the  radius  of  the  annular  cathode  rc,  i.e.  rL=rc.  From  the 
relations 

kc=N0(l0/y0,  vL=(Q0/y0)rL  and  y'0  =  (l -»i/c2)~1/2 

we  obtain 


yo  =  [l-(*rc/Wo)2]‘1/2>  n0='^  and  e o=kc 

Using  the  specific  values  1-3 < rc <  l-5cm,  a=l-98cm  and  fca  =  5-14  for  N0  =  16 
(and  ka  =  3-43  for  N0= 8),  the  operating  voltage  K0,  magnetic  field  B0,  and  cut-off 
frequencies  /0  are  determined  to  be  12<K0<16kV  (21-5 <  K0<21  kV),  0-0284  T 
(0-0387  T)  and  12-4  GHz  (8-27  GHz),  respectively.  In  the  actual  experiments  with 
the  constructed  waveguides,  signals  corresponding  to  both  2 n-  and  rc-mode 
operation  in  the  expected  frequency  regions  are  detected.  Presented  in  Figs  4  and 
5  are  spectral  outputs  recorded  by  the  HP  spectrum  analyser  (8569B).  The  signal 
at  12-5  GHz  shown  in  Fig.  4  is  the  output  of  27t-mode  radiation,  while  the  signal  at 
8-33GHz  shown  in  Fig.  5  is  the  output  of  n-mode  radiation.  In  both  cases,  the 
radiation  frequencies  are  slightly  above  the  cut-off  frequency  of  the  respective 
modes.  We  have  also  operated  the  tube  at  about  the  same  beam  voltage  and 
current,  i.e.  16  kV  and  0-5  A,  respectively.  It  is  noted  that  the  above-mentioned 
operating  voltages  are  calculated  for  an  ideal  cusp  field.  In  the  actual  experiment, 
the  cusp  field  is  far  from  ideal  and  slightly  asymmetric.  Thus,  the  16  kV  electron 
beam  can  still  pass  through  the  cusp  field  with  a  higher  downstream  magnetic  field 
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Figure  5.  Spectrum  analyser  trace  of  n-mode  radiation  at  8-33  GHz,  V  —  16kV,  and 

7=0-5  A. 


Figure  6.  Spectrum  analyser  trace  of  radiation  at  10-0 GHz,  F=  16  kV,  and  7  =  0-5  A. 


for  the  7r-mode  case  than  the  corresponding  cut-off  magnetic  field  for  the  ideal 
case.  It  is  also  noted  that  these  two  modes  occur  at  different  magnetic  fields  and 
can  be  separated  easily.  The  fact  that  only  single  line  at  each  frequency  is 
generated  suggests  that  the  waveguide  structure  can  indeed  minimize  the  compe¬ 
tition  between  the  modes.  The  spurious  signal  at  10-0 GHz  shown  in  Fig.  6  is  also 
excited  for  this  waveguide.  It  can  also  be  separated  from  27i-mode  signal  by  a 
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slight  change  in  the  downstream  guiding  magnetic  field.  At  present,  this  signal  is 
tentatively  frequency  matched  to  the  1 5th  harmonic  due  to  the  peniotron  mode. 


4.  Discussion  and  conclusion 

In  the  present  operation,  it  is  noted  that  the  detected  signal  at  12-5  GHz, 
although  distinct,  is  relatively  weak.  The  strongest  radiation  is  at  8-33  GHz,  and  at 
10  GHz  the  power  level  is  less.  It  should  be  noted  that  the  detected  signal  power 
level  does  represent  the  real  output  power.  This  is  due  mainly  to  the  inadequacy  of 
the  coupling  system  that  guides  the  radiation  from  the  circular  resonant  structure 
to  the  rectangular  receiving  antenna  horn.  In  the  present  set-up,  the  multivane 
waveguide  has  a  circular  cross-sections  of  radii  a  =1-98  cm  and  b  =  2-574  cm.  The 
structure  sits  inside  a  chamber  of  radius  5  08  cm.  Between  the  termination  of  the 
resonant  structure  and  the  end  view  port  is  a  high  voltage  break  of  5"  length 
which  is  installed  to  monitor  the  beam  current  in  order  to  determine  the  actual 
input  beam  power.  This  mismatch  in  dimension  will  cause  a  large  reflection  of 
radiated  power  back  either  through  the  waveguide  or  the  gap  between  the 
waveguide  and  the  chamber.  The  second  mismatch  is  in  the  viewport,  which  is 
0-571  cm  thick.  The  manufacturer-listed  transmission  range  for  0-32  pm  <  X  <  2-7  pm 
is  93%.  We  would  expect  the  transmission  to  be  smaller  in  the  centimeter  range  in 
which  we  currently  operate.  The  third  mismatch  is  due  to  the  use  of  an  x-band 
rectangular  horn  to  receive  the  output  signal.  In  addition  to  the  large  cross-section 
difference  in  the  radiating  aperture  and  the  receiving  horn,  this  mismatch  will 
effectively  attenuate  the  signal  prior  to  detection.  Moreover,  in  the  actual 
experiment,  when  the  horn  is  rotated  about  45°  to  the  vertical  polarization,  it  is 
discovered  that  the  detected  8-3  GHz  signal  is  maximized.  From  this  observation, 
we  can  identify  the  polarization  direction  of  this  signal.  However,  this  polarization 
is  different  for  different  modes.  This  strongly  suggested  that  there  is  a  need  to 
design  a  circular-to-rectangular  waveguide  coupler  in  order  to  optimize  the  gain  of 
the  receiver.  Therefore,  these  factors  must  be  considered  to  calibrate  the  detected 
signals  in  order  to  determine  the  actual  output  power. of  each  mode.  Nevertheless, 
after  the  calibration  the  detected  signals  are  still  lower  than  the  theoretical 
prediction.  There  are  several  candidates  which  may  contribute  to  the  discrepancy. 
The  most  dominant  factors  are: 

(1)  the  finite  thickness  of  the  annular  cathode,  which  will  cause  an  axial  energy 
spread  of  the  E-layer; 

(2)  the  non-ideal  cusp  field  will  lead  to  the  problem  of  radial  broadening  and 
coherent  off-centring  of  the  E-layer  (non-axisymmetric  E-layer). 

Another  factor  which  affects  the  overall  operating  efficiency  of  the  device  is  that  a 
large  percentage  of  the  electrons  are  lost  during  passage  through  the  anode  plate. 
This  is  an  inherent  problem  of  the  cusp  field  beam  forming  system.  In  order  to 
achieve  a  narrower  cusp  transition  width,  an  iron  plate  is  constructed  with  an 
annular  slit  which  is  slightly  larger  than  the  radial  thickness  of  the  cathode.  As  a 
result,  radial  spread  in  the  electron  Larmor  radius  may  lead  to  the  loss  of 
electrons  through  collision  with  the  anode  plate  even  before  the  beam’s  entry  into 
the  downstream  side  of  the  system.  Experimentally,  more  than  50%  of  the 
electrons  are  lost  through  the  anode. 

In  conclusion,  we  have  detected  radiation  signals  at  the  designed  frequencies  at 
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both  the  2 n  and  n  inodes.  A  signal  at  100  GHz  is  also  excited  and  has  been 
frequency  matched  to  the  15th  harmonic.  The  detected  power  level  of  the  n  mode 
is  much  higher  than  that  of  the  27:  mode,  but  the  overall  detected  power  level 
remains  lower  than  that  predicted  theoretically,  owing  mainly  to  the  poor 
receiving  unit.  In  order  to  determine  the  actual  radiated  power,  we  should  have  to 
account  for  the  losses  due  to  the  explained  factors  above.  A  comparison  of  the 
cross-section  of  the  receiving  horn  with  the  window  will  add  10 dB  to  the  detected 
signal.  An  estimated  20  dB  loss  is  attributed  to  the  mismatch  in  the  dimensions  of 
the  waveguide  and  the  chamber.  A  probable  10  dB  loss  is  assigned  to  the  reflection 
by  the  window.  Finally,  a  15  dB  loss  is  added  to  account  for  the  difference  in 
geometry  of  the  radiating  and  receiving  units.  Therefore  a  total  loss  of  55  dB  must 
be  taken  into  consideration  to  estimate  the  actual  power  output.  For  2^-mode 
operation,  the  spectrum  analyser  trace  gives  a  signal  level  of  — 18  dBm.  With  55  dB 
added,  the  estimated  radiated  power  is  37  dBm,  or  5W.  Compare  this  with  the 
theoretical  prediction  of  6-76%  for  a  16  kV  0-5  A  beam  (experimental  parameters), 
i.e.  540  W.  The  corrected  efficiency  is  still  low.  For  n-mode  operation,  the 
estimated  power  is  ~l-26kW  which  compares  favorably  with  a  theoretical 
prediction  of  2-4  kW  (30%  efficiency)  in  view  of  the  less  than  ideal  operating 
conditions  as  explained  earlier.  For  the  10  GHz  signal,  the  estimated  power  is 
about  500  W  for  the  beam  power  used  in  other  operations.  It  is  noted  that  in  the 
work  of  Chojnacki  el  al.  (1987),  the  signal  level  from  their  spectrum  analyser  trace 
is  about  the  same  as  ouf  Jt-mode  signal.  However,  Chojnacki  et  al.  (1987)  claimed 
a  power  level  of  9-1  kW.  In  order  to  verify  the  results  and  to  eliminate  the  rough 
estimation  we  have  employed  here,  a  better  coupling  system  is  required.  This  can 
be  done  by  connecting  a  coupler  having  a  cylindrical  front  end  gradually  tapered 
to  a  rectangular  waveguide  directly  onto  the  resonant  structure.  A  double-stub 
tuner  should  also  be  installed  in  the  connection  between  the  waveguide  and  the 
coupler  to  minimize  the  reflection  incurred  by  the  discontinuity  in  the  structures. 
Nevertheless,  in  the  present  operating  conditions  in  terms  of  signal  controllability 
and  step  tunability  this  device  has  shown  promise  as  a  viable  radiation  source. 
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Abstract.  Intense  magnetosonic  waves,  originally 
propagating  at  the  right  angle  with  the  interplanetary 
magnetic  field,  can  excite  a  purely  growing  mode  along 
the  interplanetary  magnetic  field  together  with  two 
symmetric  magnetosonic  sidebands  propagating  obliquely 
across  the  magnetic  field.  This  instability  process  leads 
to  the  filamentation  of  the  magnetosonic  pump  waves. 
These  two  excited  magnetosonic  sideband  modes  propagate 
together  perpendicularly  across  the  magnetic  field  and, 
meanwhile,  form  a  standing  wave  pattern  along  the 
magnetic  field.  The  thresholds  of  this  filamentation 
instability  can  be  exceeded  in  the  solar  wind  environment. 
It  is  predicted  that  the  density  fluctuations  produced  by 
the  filamentation  instability  along  the  interplanetary 
magnetic  field  have  wavelengths  greater  than,  at  least,  a 
few  Earth  radii.  The  polarization  of  the  obliquely 
propagating  magnetosonic  waves  excited  by  the 
filamentation  instability  is  determined  by  the 
characteristics  of  the  magnetosonic  pump  waves  and  the 
environmental  plasmas. 

Introduction 

Recently,  there  has  been  a  considerable  effort  in  the 
study  of  the  stability  of  Alfven  waves  in  the  space 
plasmas.  It  is  attributed  to  the  frequent  appearance,  as 
detected  by  the  space  crafts,  of  large-amplitude  Alfven 
waves  in  the  high-speed  streams  of  the  solar  wind 
[Abraham-Shranner  and  Feldman,  1977],  in  the  upstream 
Jovian  bow  shocks  [Goldstein  et  al.,  1985],  and  in  the 
interplanetary  shocks  [Vinas  et  al.,  1984],  In  general, 
these  waves  appear  to  be  circular  polarization  and 
propagating  along  the  background  magnetic  field. 
Moreover,  the  ISEE  1  and  2  spacecrafts  have  detected 
both  magnetosonic  waves  and  Alfven  waves  propagating 
upstream  into  the  solar  wind  with  phase  velocities  much 
less  than  the  solar  wind  velocity.  Hence,  these  waves  are 
blown  backward  over  the  spacecraft  and,  consequently,  an 
apparent  change  in  polarization  of  waves  occurs. 
Intrinsically  right  (left)  hand  polarized  waves  are  seen  to 
have  an  apparent  left  (right)  handed  polarization  in 
spacecraft  records  [Hoppe  and  Russell,  1983]. 

The  parametric  decay  instability  and  modulation 
instability  of  large-amplitude,  circularly  polarized  Alfven 
waves  propagating  in  the  solar  wind  along  the  magnetic 
field  have  been  analyzed  extensively  by  Lashmore-Davies 
[1976],  Goldstein  [1978],  Derby  [1978],  Wong  and  Goldstein 
[1986],  Longtin  and  Sonnerup  [1986],  and  Terasawa  et  al 
[1986],  The  filamentation  instability,  competing  with  these 
instabilities,  can  break  up  the  large-amplitude  Alfven 
waves  and  induce  not  only  plasma  density  fluctuations  but 
magnetostatic  fluctuations  [Kuo  et  al.,  1988].  Starting 
from  small  perturbations  in  plasma  density,  the 
filamentation  instability  gives  rise  to  a  modulation  of  the 
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plasma  dielectric  constant  and  then  wave  distribution  that, 
in  turn,  enhances  the  density  perturbations. 

The  modulation  instability  has  also  been  analyzed  in 
the  region  that  describes  the  nonlinear  evolution  of  the 
Alfven  wave  propagating  along  a  background  magnetic 
field.  It  is  shown  that  this  evolution  can  be  governed  by 
a  single  "derivative  nonlinear  Schrodinger  equation* 
[Mjolhus,  1976;  Spangler,  1986;  Hada  et  al,  1989],  The 
possible  applications  of  Alfven  solitons  to  solar  and 
astrophysical  plasmas  have  been  discussed  by  Ovenden  et 
al  [1983]  and  Spangler  [1985]. 

While  significant  research  has  been  directed  to  the 
Alfven  waves,  the  magnetosonic  waves  have  received  much 
less  attention.  The  magnetosonic  waves  are  hybrid  and 
elliptically  polarized  modes  that  can  propagate  obliquely 
across  the  magnetic  field.  They  become  decoupled  from 
the  Alfven  waves  when  the  direction  of  their  propagation 
is  exactly  perpendicular  to  the  ambient  magnetic  field. 
The  purpose  of  the  present  work  is  to  investigate  the 
filamentation  instability  of  the  magnetosonic  waves  that 
propagate  perpendicularly  across  the  magnetic  field.  This 
instability  excited  obliquely  propagating  magnetosonic 
sideband  modes.  The  concomitantly  excited  nonoscillatory 
(purely  growing)  modes  have  their  wave  normal  parallel  to 
the  magnetic  field.  By  contrast,  the  purely  growing 
modes  excited  by  the  filamentation  instability  of  Alfven 
waves  are  magnetic  field-aligned  [Kuo  et  al.,  1988,  1989]. 
Hence,  the  characteristics  of  the  filamentation  instability 
under  consideration  are  distinctively  different  from  those 
of  the  instability  discussed  in  our  previous  work  [Kuo  et 
al.,  1988,  1989]. 

This  paper  is  organized  as  follows.  The  coupled 
mode  equations  are  derived  in  section  2.  They  show  how 
the  magnetosonic  sidebands  and  the  purely  growing  modes 
are  parametrically  coupled  through  the  magnetosonic  pump 
wave.  In  section  3  the  dispersion  relation  is  obtained  and 
analyzed  to  derive  the  threshold  conditions  and  the 
growth  rates  of  the  filamentation  instability.  The 
relevance  of  the  present  work  to  some  observations  is 
discussed.  Finally  presented  in  section  4  are  summary  and 
conclusions. 

Filamentation  Instability  of  Magnetosonic  Waves 

Consider  a  magnetosonic  wave  propagating 
perpendicularly  across  the  interplanetary  magnetic  field 

B0=zB0.  Its  wave  magnetic  field  is  represented  by 

Bp  -  zBpexplKkoX-Wot)]  +  c.c.  ( 1) 

where  the  unperturbed  wave  field  intensity  Bp  is 
assumed  to  be  constant  and  a  real  quantity  for  simplicity. 
The  wave  number  kg  and  the  angular  frequency  w0 
are  related  by  the  dispersion  relation 
A+c*t),  wherein  vA  and  c,  are  the  Alfven  speed 
and  the  ion  acoustic  speed,  respectively.  The  velocity 
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responses  of  electrons  and  ions  to  the  magnetosonic  pump 
wave  can  be  written  as 

Vp«  =  [x+i(kJ0v2A/o)0q)y](oJ0/k0)(Bp/B0)exp(i(k0x-w0t)] 

+  c.c. 

Vpi  =  x(  u0/  !<„)  ( Bp/B0)  exp(  i(  k<,x-w0t)  ]  +  c.c.  (2) 

In  the  following  we  will  investigate  a  filamentation 
instability  which  can  break  up  a  large-amplitude 
magnetosonic  wave  into  filaments.  The  magnetosonic  wave 
is  assumed  to  be  initially  uniform  as  described  by  (1)  for 
simplicity.  The  filamentation  instability  results  from  small 
fluctuations  in  plasma  density  which  lead  to  modulation  of 
plasma  dielectric  constant  and  spatial  wave  distribution, 
which,  in  turn,  enhance  the  plasma  density  fluctuations. 
Such  a  positive  feedback  process  involves  a  four- wave 
coupling  process:  through  the  magnetosonic  pump  wave 

(S’ko)>  two  side-band  perturbations  (w+)  are 
coupled  with  a  purely  growing  density  perturbation 
(w,kz).  The  wave  frequency  and  wave  vector 
matching  conditions,  are  u>+-of=u0=w_+o/  and  1c  +="k0±"k,  are 
satisfied  where  the  asterisk  denotes  the  complex 
conjugate.  Illustrated  in  Figure  1  is  the  wavevector 
matching  condition,  showing  the  relative  orientation  of 
these  wave  vectors.  It  will  be  sown  that  these 
perturbations  can  be  excited  simultaneously.  In  other 
words,  they  grow  exponentially  with  the  same  growth 
rate,  7=-iw  at  the  expense  of  magnetosonic  pump  wave 
energy.  This  instability  requires  a  threshold  to  occur, 
which  is  going  to  be  determined  later.  In  the  analysis  of 
filamentation  instability,  pump  field  intensity  is  assumed 
to  be  a  constant  maintained  by  a  source.  Pump  depletion 
effect  which  may  contribute  to  the  saturation  of 
filamentation  instability  is,  therefore,ignored  within  the 
present  framework  of  instability  analysis. 

Our  proposed  theory  is  developed  on  the  basis  of  the 
two-fluid  plasma  model.  The  fluid  model  is,  in  general, 
valid  provided  that  the  condition  k2<^/f^«l  is 
satisfied  where  k,  c,,  and  are  the  wave  number 
of  the  purely  growing  mode,  the  ion  acoustic  speed,  and 
the  ion  gyrofrequency,  respectively.  Therefore,  the 
analysis  of  filamentation  instability  is  limited  to  the 
region  of  k/  ^<1.2.  To  analyze  the  instability. 


A 
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Fig.  1.  The  relative  orientation  and  wave  vector  matching 
relation  for  the  magnetosonic  pump  (coo.ko),  two  sidebands 
(co±)k±)  and  the  purely  growing  mode  (co,k)  as  required  in  the 
proposed  four-wave  coupling  process. 


the  complete  set  of  fluid  equations  will  be  linearized  with 
respect  to  the  first-order  perturbations  in  physical 
quantities  such  as  density  (n),  velocity  (v),  wave  electric 

field  (?)  and  wave  magnetic  field  (?).  These 
first-order  physical  quantities  associated  with  the  sideband 

modes  are  denoted  by  6n±,  6"v  t^±,6?±,  and  SB±  those 
associated  with  the  purely  growing  mode  n,, 

V...,  ?t,  and  Thus,  the  total  density  is  given 

by  n=n0+n,  +  (  Sn+  +  fn_+c.  c. ) ;  similarly, 

*(  .4+5V.4++i"%i-+c-c->  Ep+«E++«E.+c.c. ) , 

and  B=B0+(  Bp+5B++SB_+c.c.)  for  the  concerned  four-wave 
interaction  process. 

The  purely  growing  mode  is  associated  with 
electrostatic  disturbances,  varying  spatially  along  the 

magnetic  field,  i.e.,  k=zk.  The  fluid  equations 
used  to  derive  the  coupled  mode  equation  for  the  purely 
growing  mode  include  the  linearized  continuity  equations 
for  electrons  and  ions: 

If  +  “o  |  v.„  -  0  -  |  n.+no  |  vgil  (3) 

and  the  combined  linearized  equations  of  motion  for 
electrons  and  ions: 

m.  |  v.«  +  mi  |  v.i.  +  <F« 

+  F„)  -  -  m,  (cj/n,,)  |  n. 

In  deriving  the  above  equations,  quasi -neutrality  has  been 
assumed.  Ftl  and  Fif  represent  the  z  components  of  the 
nonlinear  Lorentz  forces  experienced  by  electrons  and 
ions,  respectively,  namely,  F,jt«<-mtj(  •  Vvtp 
+  q^Cvej/c)  xB>,  where  angle  brackets  denote  the 
time  average  over  the  pump  wave  period.  The  explicit 
expressions  of  F,^,  are  given  by 

=  |  •  ^.i> 

+i£«i  [*k)  I  ‘  <4> 

where  S~v t^S~v is  the  velocity  responses  of 
electrons  and  ions  to  the  combined  sideband  mode; 
^(^)=+ !(-!)•  The  expression,  (4)  reduces  to  the 
form  of  ponderomotive  force  in  an  unmagnetized  plasma. 
Without  losing  the  generality,  the  density  fluctuations  of 
the  purely  growing  mode  can  be  expressed  as 

n,  -  n,  exp(Tt)cos  kz 


where  n,  and  7  are  the  real  amplitude  and  the 
growth  rate  of  the  density  fluctuations,  respectively. 

Combining  the  equations  in  (3)  yields  the  coupled 
mode  equation  for  the  purely  growing  mode  as 

(7J+k1c2<)(n1/n0)  -  (k/mi)(Pn+Fit)  (5) 
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where  F.^»F,^exp(it)  sinkz  is  assumed.  One  can 
see  from  (4)  and  (5)  that  the  density  fluctuations 
associated  with  the  purely  growing  mode  stem  from  the 
differential  interaction  of  the  magnetosonic  pump  wave 
and  sideband  modes.  Hence,  expressions  for  the  physical 
quantities  of  the  combined  sideband  mode  have  to  be 
determined  in  a  self-consistent  manner.  The  purely 
growing  mode  is  realized  to  be  a  nonlinearly  driven  mode. 
However,  in  the  absence  of  the  nonlinear  Lorentz  force, 
(5)  reduces  to  the  linear  dispersion  relation  of  the  ion 
acoustic  wave. 

The  magnetosonic  sideband  modes  are  excited  through 
the  beating  current  driven  by  the  magnetosonic  pump 
wave  field  in  the  density  fluctuations  of  the  purely 
growing  mode.  It  should  be  noted  that,  since  the 
magnetosonic  pump  wave  propagates  perpendicularly  across 
the  magnetic  field  and  the  purely  growing  mode  oscillates 
spatially  along  the  magnetic  field,  two  magnetosonic 
sidebands  are  excited  symmetrically  around  the  pump  wave 
and  they  propagate  obliquely  across  the  magnetic  field 
(see  Figure  1).  Thus,  these  two  sideband  modes  propagate 
together  in  the  propagation  direction  of  the  pump  wave 
and  form  a  standing  wave  pattern  along  the  magnetic 
field.  In  the  following  derivation  of  the  coupled  mode 
equations  for  the  two  sideband  modes,  we  can,  therefore, 
combine  them  into  a  single  coupled  mode  equation  for  the 
combined  sideband  mode. 

This  combined  sideband  mode  is  characterized  by  the 
following  physical  quantities  with  self-consistent 
expressions: 

Sn  =  Sn  exp(Tt)cos  kz  exp(i( l^x-Uot) ]  +  c.c. 

£vtj  m  +  iy«v#iy)cos  kz  +  iz«v,^  sin  kz]exp(Kt) 

•  expfi(kgx-w0t)]  +  c.c. 

$1?  -  [(k/koH-ix+ayJsin  kz  +  z  cos  kz]5B  exp(Tt) 

•  exp[i(  k„x-<D0t)  ]  +  c.c.  (6) 

where  a,  an  unknown  parameter  defining  the  polarizations 
of  the  sideband  modes,  will  be  determined  later  on.  The 
fluid  equations,  used  to  derive  the  coupled  mode  equation 
for  the  combined  sideband  mode,  include  the  linearized 
continuity  equations  for  electrons  and  ions: 

|  («n/n0)  +  V  •  [fiV,  +  (n,/n0)'vp<]  - 
0-f  (Sn/n,,)  +  V  .  +  (n,/n0)'vpi]  (7) 

the  linearized  Maxwell  equation  for  induction  law: 

|  («B/B0)  -  V  x  (ff?,  x  z)  (8) 

the  combined  equations  of  motion  for  electrons  and 
ions: 

5  .  -  c3,  V  ( fin/n 0)  +  r\(rvrf?,)  X  z  (9) 

and  the  linearized  Maxwell  equation  for  Amperes'  law: 

V  x  («1/B0)  -  (tyv^K^-fV.) 

+  (n,/n0)(Vpi  -  Vp.)]  (10) 


Substituting  (6)  into  equations  (7) -(10),  we  obtain 

=  (VV(®/Bo) 

«*.y  *  CtK/koMffi/Bo) 

-  -  (V^)U^0cVrfi)(®/»o)- 

(l-cw0/q)(«g/B0)]/(l-wJ0/rfi) 

^iy  -  ( Vf\M(kocV(,'oH®/no>  +  wq/H)-i] 

-  (kcJ|/oJ0)(65f/n0) 

+  o<kv2A/f\)(«g/B0) 

«/"o  “  {[(  l-«Vty/(  l-«fyrfi)K«B/B0) 

+  (K./n0)(Bp/B0))/[l-k2c2,/tv20  +  k20cV(tfru’20)]  (U) 

Using  the  results  in  (II)  together  with  the  x  and  y 
components  of  ( 10) ,  we  can  derive  the  following  coupled 
mode  equation  for  the  combined  sideband  mode: 

[aV^+^oHvVc2.)  +  k*(  k^o)  vVV^O 
+  o<w0/^)[k20c,.-wJo(  l-k2c2,/^0)])(«B/B0) 

-  -  l-^./^o)  -  ^.KVnoMBp/Bo)  (12) 

as  well  as  the  expression  for  the  parameter  ac 

a  -  (VtVll-^k,+ka0)cV<*’20]/ 

[  l-cs,/vsA-(  kMA/o*0)  ( l-kV./^o)  ]  ( 13) 

In  the  derivation  of  (12)  and  (13),  <*’*0/n,i  « 
1  has  been  assumed.  Equation  ( 12)  shows  that  the 
magnetosonic  sideband  fields  are  produced  by  the 
magnetosonic  pump  field  in  the  density  fluctuations  of  the 
purely  growing  mode.  Equations  (5)  and  (12)  thus  form 
the  complete  set  of  coupled  mode  equations,  which  can  be 
based  on  to  analyze  the  proposed  filamentation  instability 
of  magnetosonic  waves  in  the  solar  wind  environment. 

Analysis  of  the  Proposed  Instability 

The  dispersion  relation  of  the  concerned  instability  can 
be  derived  from  the  two  coupled  mode  equations,  (5)  and 
(12).  Its  derivation  is  briefly  described  as  follows.  With 
the  aid  of  (2),  (4),  (6),  and  (11),  the  coupling  term  on  the 
right-hand  side  of  (S)  can  be  expressed  as 

k(  F.+PjJ/mj  -  l-aaj0/q)/ 

( l-k*c*,/w*| 0J(Bp/B0)(4§/B0  +  «g‘/B0) 

-  2[kJc3,/(  l-l^c* ,/<»20)](ff</n0)( Bp/Bo) *  (14) 

Substituting  (14)  into  (5)  and  combining  it  with  (12)  yield 
the  dispersion  relation 
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Fig.  2.  a/8  versus  k/ko  for  several  selected  values  of  P  in  the  range  of  0.25-2.5;  it  shows  how  the 

polarization  of  the  excited  magnetosonic  sidebands  is  determined  by  the  characteristics  of  the  magnetosonic 
pump  waves  and  the  environmental  plasmas. 


T'+kV,  -  2kVA<  1 1-K  1+0  /( 0- 1 
[I-«  l+mi-oQ/^l-^/k2,,)] 

( 0  l-flc2/  k*0)  /[( 0 1)  l+k1/^)  +  a«0i)(kVk,-6] 

-( k2/  k*0)  A  01)  /( 0  l-/*2/  k2,,) }  ( Bp/B,,) 2  (15) 

where  0c2,/  v*A  is  the  ratio  of  the  plasma 
kinetic  pressure  to  the  background  magnetic  pressure,  and 
£=“0/ l^- 

Setting  7=0  in  ( 15)  determines  the  threshold  field  of 
the  filamentation  instability 

(Bpth/B0)1  -  (k*/2l^)/t(01)*-A  l+k2/^) 

+  ofl(*i)(kVi^A]/{M*i)/ 

(01-A*/k*o][ 1-fl 0- 1) ( l-crf) /( 0  1-/5IC2/ k2,) ] 

OM-ZSk2/^)-  (k*/k?o)A01){(/M)s  -  fll+k2/^) 

+  a«(01)(k2o/k2-0]/(01-/k2/k2o)}  (16) 

In  terms  of  (16),  the  growth  rate  of  the  instability 
derived  from  (15)  is  found  to  be 

i  -  [( VBPt»>J-1l1/,kc.  (,7> 

The  threshold  condition  defined  by  (16)  is 
conveniently  expressed  for  the  normalized  threshold  field 
intensity  (Bpth/B0)  as  a  function  of  the 

normalized  wave  number  (k/kg).  The  parameters 

P{  -c2,/  v*A) ,  5(=<^,/I^),  and  a  (defined 

by  (13))  characterize  the  solar  wind  plasmas  and  the 


propagating  magnetosonic  waves.  In  fact,  the  parameter 
a  can  also  be  expressed  in  terms  of  (k/kg),  5,  and  P 
as  follows: 

a  -  5n-2fl+(k/kg)V/(  1+0/1  l-0[(k/kg)2/ 

( 1+0  ][l-(  k/kg)  20(  1+01)  (18) 

Displayed  in  Figure  2  is  (a/0  as  a  function  of  (k/kg) 
for  several  values  of  P  chosen  from  the  range  of  0.25- 
2.5.  The  result  shows  that  a  is,  in  general,  larger  than 
one  for  5<0.3  which  is  the  region  of  practical 
interest.  Especially,  a»l  for  01  and  5«1.  In  this 
parameter  range,  the  decay  sidebands  become  linear 
polarization  on  the  y-z  plane  and  are  dominated  in  the 
nature  of  shear  mode. 

The  threshold  fields  (Bptb/B0)  as  a  function 
of  (k/kg)  are  shown  in  Figures  3-6.  For  illustrative 
purposes,  the  values  of  P  are  again  selected  from  the 
range  of  0.25-2.5  and  5=0.01,  0.03,  0.06,  or  0.3  is  chosen. 
The  corresponding  growth  rates,  for  ( Bp/B0)2=0.05, 
as  a  function  of  k  are  given  in  Figures  7-10.  Several 
outstanding  features  of  the  proposed  instability  can  be 
seen  from  these  figures.  For  instance,  the  threshold  field 
generally  increases  with  k,  the  wave  number  of  the 
excited  purely  growing  mode,  and  it  finally  reaches 
infinity  at  certain  k<kg.  While  the  threshold  field 
increases  with  P  in  the  region  O<01,  it  decreases  with 
P  for  P>  1.  The  growth  rate  has  a  peak  value  (denoted 
by  <1  )  at  certain  k  (denoted  by  kp,).  We 

note  that  both  7m|o(  and  kp,  decrease  with  P 
for  O<01,  but  they  increase  with  P  for  P>  1.  By 

contrast,  while  7 decreases  with  5,  k_ 

increases  with  5.  The  thresholds  and  growth  rates  of 
the  instability  vary  with  the  magnetosonic  wave  frequency. 
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Fig.  3.  The  normalized  threshold  field  intensity  (Bp/Bo)  versus  (k/ko)  with  several  selected  values  of  P  in  the 
range  of  0.25-2.5  for  8=0.01. 


It  is  found  that  the  threshold  fields  (growth  rates) 
increase  (decrease)  with  the  magnetosonic  wave 
frequency. 

Let  us  base  upon  the  observed  characteristics  of 
magnetosonic  waves  (Hoppe  and  Russell,  1983,  and 


references  therein]  to  examine  the  excitation  of  the 
proposed  filamentation  instability.  The  magnetosonic 
waves  observed  by  the  ISEE  1  and  2  typically  have  low 
frequencies  of  order  of  0.01  Hz,  which  is  about  0.1  times 
the  local  proton  gyrofrequency  (i.e. ,  5-uo/f}-0. 1) . 


Fig.  4.  (Bp/Bo)  versus  (k/ko)  for  8=0.03. 
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Fig.  5.  (Bp/Bo)  versus  (k/ko)  for  8=0.06. 


The  local  interplanetary  magnetic  field  (B0),  therefore, 
is  of  order  of  100  nT.  The  magnetosonic  waves  typically 
have  amplitudes  ( Bp)  of  about  10  nT,  namely, 
Bp/B0  ~  0. 1  with  wavelengths  ~1  RE.  The  parameters 
adopted  in  the  analysis  of  our  proposed  instability  are 
generally  in  agreement  with  these  observations.  Under 


the  condition  of  Bp/B0  ~  0.1,  the  thresholds  of 
the  proposed  instability  can  be  exceeded  for  k/ kg 
(sAg/A)  less  than  0.1  or  so,  where  \  and  A  are  the 
wavelength  of  the  magnetosonic  wave  and  the  scale  length 
of  the  excited  purely  growing  mode,  respectively.  An 
important  prediction  of  the  proposed  process  is  that 


Fig.  6.  (Bp/Bo)  versus  (k/ko)  for  8=0.3. 
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plasma  density  fluctuations  with  scale  lengths  longer  than 
at  least,  a  few  RE  (Earth  radius)  can  be  generated  along 
the  magnetic  field  by  magnetosonic  waves.  The 
polarization  of  the  obliquely  propagating  magnetosonic 
waves  is  determined  by  the  value  of  or  which  in  turn,  is 
a  function  of  the  characteristic  parameters  of  the 
environment:  S,  P,  and  kg  (see  equation  (18)). 


However,  we  note  that  there  seems  to  exist  a  discrepancy 
between  the  observed  and  the  predicted  propagation 
directions  of  magnetosonic  waves.  The  waves  upstream  of 
Earth's  bow  shock  with  frequencies  of  0.01  Hz  were 
Observed  to  propagate  nearly  parallel  to  the  magnetic 
field,  while  the  predicted  magnetosonic  waves  propagate 
obliquely  across  the  field.  We  suspect  that  this 


Fig.  8.  y  versus  (k/ko)  for  8=0.03. 
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Fig.  9.  Y  versus  (k/ko)  f°r  8=0.06. 


discrepancy  may  arise  from  the  intrinsic  difficulty  in  the 
in  situ  measurements  of  upstream  waves  by  satellites. 
Since  the  speed  of  solar  wind  is  far  greater  than  those  of 
satellites,  all  measures  waves  appear  to  be  nearly  parallel 
propagating.  The  predicted  oblique  propagation  of 
magnetosonic  waves  excited  by  filamentation  instability 
awaits  the  corroboration  of  future  observations. 


4.  Summary  and  Conclusions 

The  excitation  of  a  parametric  instability  by 
magnetosonic  waves  in  the  solar  wind  environment  is 
proposed.  The  magnetosonic  pump  wave  originally 
propagates  at  the  right  angle  with  the  magnetic  field. 
This  instability  can  excite  a  purely  growing  mode  along 


Fig.  10.  Y  versus  (k/ko)  for  8=0.3. 
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the  magnetic  fieid  together  with  two  symmetric 
magnetosonic  sidebands  propagating  obliquely  across  the 
magnetic  field.  These  two  sideband  modes  thus  propagate 
together  perpendicularly  across  the  magnetic  field  and, 
meanwhile,  form  a  standing  wave  pattern  along  the 
magnetic  field.  The  physical  process  responsible  for  the 
instability  can  be  described  as  the  scattering  of  the 
magnetosonic  pump  wave  into  obliquely  propagating 
sidebands  by  the  simultaneously  excited  plasma  density 
fluctuations,  wnich  are  associated  with  a  purely  growing 
mode.  This  instability  leads  to  the  filamentation  of  the 
magnetosonic  pump  wave.  The  thresholds  of  this 
filamentation  instability  can,  indeed,  be  exceeded  in  the 
solar  wind  environment.  The  density  fluctuations 
produced  by  the  filamentation  instability  along  the 
magnetic  field  are  predicted  to  have  wavelengths  greater 
than,  at  least,  a  few  Earth  radii.  As  shown  in  (18),  the 
polarization  of  the  obliquely  propagating  magnetosonic 
waves  excited  by  the  filamentation  instability  is 
determined  by  the  environmental  characteristic  parameters: 
«(-<*b/H),  0(  “C1*/  v*A)  ,  and  k„.  The 

predicted  oblique  propagation  of  magnetosonic  waves, 
however,  has  not  been  confirmed  by  the  available 
observations  yet. 
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Bragg  scattering  of  electromagnetic  waves  by  microwave-produced 
plasma  layers 

S.  P.  Kuo  and  Y.  S.  Zhang 

Weber  Research  Institute,  Polytechnic  University,  Route  110,  Farmingdale,  New  York  11735 
(Received  18  Juiy  1989;  accepted  26  October  1989) 

A  set  of  parallel  plasma  layers  is  generated  by  two  intersecting  microwave  pulses  in  a  chamber 
containing  dry  air  at  a  pressure  comparable  to  the  upper  atmosphere.  The  dependencies  of 
breakdown  conditions  on  the  pressure  and  pulse  length  are  examined.  The  results  are  shown  to 
be  consistent  with  the  appearance  of  tail  erosion  of  the  microwave  pulse  caused  by  air 
breakdown.  A  Bragg  scattering  experiment,  using  the  plasma  layers  as  a  Bragg  reflector,  is 
then  performed.  Both  time  domain  and  frequency  domain  measurements  of  wave  scattering  are 
conducted.  The  experimental  results  are  found  to  agree  very  well  with  the  theory. 


I.  INTRODUCTION 

Over  the  years  there  has  been  considerable  effort  con¬ 
centrated  on  investigating  the  propagation  of  high-power 
microwave  pulses  through  the  air.  Ionization  enhances  the 
electron  density  in  the  plasma  trail  following  the  pulse  and  in 
turn,  radically  modifies  wave  propagation.  ,"6  However,  the 
possibility  of  generating  plasma  layers  by  high-power  radio 
waves  in  the  atmosphere  at  heights  of  30-60  km  was  first 
studied  theoretically  by  Gurevich  and  his  co-workers.7"10 
An  artificial  atmospheric  plasma  is  thought  to  be  capable  of 
reflecting  radio  waves  and  extends  the  capability  of  iono¬ 
spheric  radio  communications.  In  general,  there  are  two  pro¬ 
cesses  that  can  be  used  for  the  reflection  of  radio  waves.  One 
process  is  through  the  cutoff  reflection  for  which  the  fre¬ 
quency  of  the  incident  wave  has  to  be  lower  than  the  maxi¬ 
mum  effective  plasma  frequency.  Therefore  the  maximum 
frequency  of  a  radio  wave  reflecting  from  the  plasma  is  deter¬ 
mined  by  the  electron  density  of  the  plasma  and  the  collision 
frequency  between  electrons  and  neutrals.  On  the  other 
hand,  if  the  plasma  is  structured  with  layers,  it  is  believed 
that  Bragg  scattering  would  be  an  effective  process  for  re¬ 
flecting  radio  waves  off  the  plasma  layers.  In  this  case,  there 
is  no  upper  bound  in  the  frequency  of  reflection;  however, 
the  Bragg  condition  determines  the  region  of  frequency  of 
radio  waves  for  achieving  the  maximum  reflection.  One  of 
the  apparent  advantages  of  the  Bragg  scattering  process  is 
that  much  less  dense  plasma  than  that  required  for  the  cutoff 
reflection  process  can  be  used  for  reflecting  waves  of  the 
same  frequencies.  This  implies  that  much  less  power  is  need¬ 
ed  for  plasma  generation  and  its  sustainment  if  the  Bragg 
scattering  process  is  adopted.  Although  the  use  of  the  Bragg 
scattering  process  requires  the  scattering  plasma  to  be  struc¬ 
tured,  it  may  well  be  the  inherent  property  of  the  plasma 
depending  on  the  scheme  of  plasma  generation. 

It  turns  out  that  the  most  promising  method  of  plasma 
generation  by  electromagnetic  radiation  is  to  strike  a  dis¬ 
charge  in  the  intersection  region  of  parallel  polarized  wave 
beams."  The  interference  between  the  fields  of  the  two 
beams  (for  example)  enhances  the  peak  field  amplitude  and 
thus  reduces  the  required  power  level  of  each  yf  beam.  This, 
in  turn,  helps  to  reduce  the  propagation  loss  in  pulse  energy 
before  the  two  beams  intersect.  In  fact,  more  energy  will  be 


delivered  to  the  destination  because  the  pulse  tail  erosion 
problem  can  be  almost  completely  suppressed,  especially 
when  the  intersection  altitude  is  chosen  to  be  near  50  km  ( 1 
Torr  pressure)  where  the  breakdown  threshold  (for  yf 
pulses  in  the  GHz  frequency  range)  is  minimum  (so  that  the 
most  effective  ionization  can  be  achieved).  In  the  intersec¬ 
tion  region,  the  field  amplitude  varies  periodically  in  space  in 
the  direction  perpendicular  to  the  plane  bisecting  the  two 
beams.  Consequently,  the  induced  ionization  patch  consists 
of  a  set  of  parallel  plasma  layers,  as  an  ideal  set  up  for  Bragg 
scattering  purposes. 

In  this  paper,  we  report  some  laboratory  experiments  on 
the  characteristic  features  of  air  breakdown  produced  by  two 
microwave  beams  intersecting  at  right  angles  to  each  other 
and  Bragg  scattering  of  test  waves  off  the  induced  plasma 
layers.  The  organization  of  this  paper  is  as  follows.  The  ex¬ 
perimental  setup  is  described  in  Sec.  II,  in  which  the  experi¬ 
mental  results  characterizing  the  dynamics  and  structure  of 
induced  plasma  are  also  given.  In  Sec.  III  the  Bragg  scatter¬ 
ing  experiment  is  reported.  Both  time  domain  and  frequency 
domain  measurements  of  scattering  signals  are  performed. 
The  experimental  results  are  compared  with  the  theory.  The 
work  is  summarized  and  discussed  in  Sec.  IV. 

II.  CHARACTERISTICS  OF  AIR  BREAKDOWN  BY 
MICROWAVE  PULSES 

Chamber  experiments  to  study  the  Bragg  scattering  off 
microwave-generated  plasma  layers  are  conducted.  The 
chamber  is  a  2  ft  Plexiglas  cube  that  is  filled  with  dry  air  to  a 
pressure  corresponding  to  the  simulated  altitude.  The  mi¬ 
crowave  power  is  generated  by  a  single  magnetron  tube 
(OKH  1448)  driven  by  a  soft  tube  modulator.  The  magne¬ 
tron  delivers  1  MW  peak  output  power  at  a  center  frequency 
of  3.3  GHz.  The  modulator  uses  a  pulse  forming  network 
having  a  pulse  width  that  can  be  varied  from  1. 1  to  3.3  /xsec, 
with  respective  repetition  rates  from  60  to  20  Hz.  The  micro- 
wave  beams  are  fed  into  the  cube,  with  parallel  polarization 
direction,  by  two  S-band  microwave  horns  placed  at  right 
angles  to  the  adjacent  sides.  The  plasma  layers  are  then  gen¬ 
erated  in  the  central  region  of  the  chamber  where  the  two 
beams  intersect.  Shown  in  Fig.  1  is  a  photo  of  the  plasma 
layers  that  are  manifested  by  the  enhancement  of  the  airglow 
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FIG.  1.  A  top  view  of  plasma  layers  produced  in  the  chamber 


from  the  corresponding  locations.  A  maximum  of  eight  lay¬ 
ers  can  be  generated,  though  only  a  few  of  them  are  shown  in 
the  photo.  Shown  in  Fig.  2(a)  is  the  typical  envelope  of  a  1.1 
(j. sec  pulse  used  for  plasma  generation.  Using  a  focusing  lens 
to  localize  the  enhanced  airglow,  its  temporal  evolution  be¬ 
tween  the  two  consecutive  pulses  is  then  recorded  on  the 
oscilloscope  through  a  photomultiplier  tube.  A  typical  result 
is  shown  in  Fig.  2(b),  which  shows  the  growth  and  decay  of 
the  enhanced  ionization.  We  have  also  measured  the  break¬ 
down  threshold  field  as  a  function  of  the  pressure.  The  mi¬ 
crowave  field  is  measured  by  a  microwave  probe  that  has 
been  calibrated  by  a  known  waveguide  field.  The  breakdown 


i  ?jrV 


A 

v  ; 

\ 


(a) 


(b) 


FIG.  2.  The  optical  measurement  of  the  growth  and  decay  of  the  airglow 
from  a  plasma  generated  by  two  intersecting  1 . 1  ftsec  microwave  pulses.  ( a ) 
Envelope  of  a  1 . 1  fjse c  microwave  pulse;  the  horizontal  axis  is  200  nsec/divi¬ 
sion.  (b)  Growth  and  decay  of  plasma  glow  measured  by  an  optical  probe,  1 
/isec/division. 


threshold  is  defined  as  the  wave  field  that  is  able  to  enhance 
minimum  observable  airglow  monitored  by  the  photomulti¬ 
plier  tube.  Shown  in  Fig.  3  are  the  Paschen  breakdown 
curves  for  the  cases  of  1 . 1  and  3.3  /zsec  pulses.  Since  a  shorter 
pulse  requires  a  larger  ionization  rate  in  order  to  generate  the 
same  amount,  of  electrons  (which  are  proportional  to  the 
enhanced  airglow  flux),  the  threshold  field  is,  therefore,  ac¬ 
cordingly  increased.  This  tendency  is  clearly  demonstrated 
in  Fig.  3.  It  shows  that  the  breakdown  threshold  field  for  a 
1.1  /zsec  pulse  is  always  larger  than  that  for  a  3.3  /zsec  pulse. 
The  results  also  show  that  in  both  cases,  the  breakdown 
threshold  field  decreases  with  a  decrease  in  air  pressure  and 
reaches  a  minimum  in  the  2  to  1  Torr  region,  where  <u  =~  vc ;  <u 
and  vc  are  the  microwave  frequency  and  electron-neutral 
collision  frequency,  respectively.  With  a  further  decrease  in 
the  pressure,  the  breakdown  threshold  field  increases  again. 
The  increase  of  the  threshold  field  also  occurs  because  the 
wave  is  in  the  pulse  mode.  The  ionization  frequency  and 
collision  frequency  are  proportional  to  the  neutral  density; 
lower  pressure  requires  a  larger  field  in  order  to  maintain  the 
ionization  frequency. 

The  dependence  of  the  breakdown  threshold  field  on  the 
pressure  is  also  manifested  by  a  similar  dependence  of  the 
degree  of  attenuation  in  the  tail  portion  of  a  single  transmit¬ 
ted  pulse  through  the  chamber.  The  experiment  is  performed 
by  reducing  the  chamber  pressure  consecutively  from  8  Torr 
to  50  m  Torr,  while  the  incident  pulse  is  fixed  at  constant 
amplitude.  A  series  of  snapshots  demonstrating  this  behav¬ 
ior  is  presented  in  Fig.  4.  In  the  high  pressure  region  (>8 
Torr),  the  breakdown  threshold  field  is  higher  than  that  of 
the  incident  pulse  and,  therefore,  very  little  ionization  can 
occur;  thus  the  pulse  can  pass  through  the  chamber  almost 
without  any  distortion  [Fig.  4(a)].  However,  as  the  pres¬ 
sure  drops,  the  breakdown  threshold  also  decreases  before 
reaching  the  minimum,  and  hence,  more  ionization  occurs 
and  so  does  more  distortion  to  the  pulse  [Figs.  4(b)-4(d)  ] . 
The  distortion  always  starts  from  the  tail  portion  of  the  pulse 


FIG.  3.  Dependence  of  the  air  breakdown  threshold  field  on  the  pressure  for 
two  pulse  lengths. 
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FIG.  4.  Pulse  received  after  passing  through  the  chamber  at  different  air  pressures,  (a)  p  =  8Torr,  (b )  p  =  6Torr,  (c)  p  =  5Torr,  (d)/>  =  4Torr,  (e)  p  =  2 
Torr,  (0  p  =  1  Torr,  (g)  p  —  0.45  Torr,  (h)  p  =  0.2  Torr,  (i)  p  =  0.1  Torr,  (j)  p  =  0.07  Torr,  and  (k)  p  =  0.05  Torr. 


(i.e.,  tail  erosion )  because  it  takes  a  finite  time  for  the  plasma 
to  build  up  and  thus  maximum  absorption  of  pulse  energy  by 
the  generated  electrons  always  appears  in  the  tail  of  the 
pulse.  Consequently,  the  leading  edge  of  the  pulse  is  usually 
not  affected.  Between  2  to  1  Torr,  the  pulse  appears  to  suffer 
maximum  tail  erosion  and  hence  only  the  very  narrow  lead¬ 


ing  edge  of  the  pulse  can  pass  through  the  chamber  [Figs. 
4(e)  and  4(0  ].  The  tail  erosion  becomes  weak  again  for  a 
further  decrease  in  the  pressure  [Figs.  4(g)-4(j)j  and  even¬ 
tually  vanishes  [Fig.  4( k )  ]  once  the  pressure  becomes  so 
low  (<0.05)  Torr  that  the  breakdown  threshold  power  ex¬ 
ceeds  the  peak  power  of  the  incident  pulse. 
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III.  BRAGG  SCATTERING 

The  plasma  layers  generated  by  the  two  crossed  micro- 
wave  pulses  are  used  for  the  Bragg  scattering  study.  The 
spatial  distribution  of  the  plasma  layers  is  first  measured 
with  a  Langmuir  double  probe.  This  is  done  by  using  a  mi¬ 
crowave  phase  shifter  to  move  the  plasma  layers  across  the 
probe.  The  peak  density  distribution  for  a  spatial  period  is 
thus  obtained  and  presented  in  Fig.  5.  The  result  shows  that 
we  have  indeed  produced  very  sharp  plasma  layers  with  very 
good  spatial  periodicity  (Fig.  1 ).  A  Bragg  scattering  experi¬ 
ment  has  then  been  conducted  and  described  as  follows. 

Presented  in  Fig.  6  is  a  block  diagram  of  the  experiment 
setup.  In  addition  to  the  facility  used  for  plasma  generation 
( located  to  the  left  of  the  Plexiglas  chamber ) ,  a  sweep  micro- 
wave  generator  (4-8  GHz)  is  used  to  generate  a  test  wave 
that  is  incident  into  the  chamber  through  a  C-band  horn. 
The  incident  angle  of  the  test  wave  with  respect  to  the  nor¬ 
mal  of  the  plasma  layers  is  45*.  Hence  the  S-band  horn  #2 
located  at  a  right  angle  to  the  adjacent  side  can  be  used  as  the 
receiver  of  the  Bragg  scattering  test  wave.  In  order  to  sepa¬ 
rate  the  Bragg  coherent  reflection  mechanism  from  the  cut¬ 
off  reflection  mechanism,  the  test  wave  is  swept  in  a  frequen¬ 
cy  range  much  higher  than  the  plasma  cutoff  frequency. 
Consequently,  the  test  wave  will  be  received  by  the  S-band 
horn  #  1  even  while  the  plasma  is  present.  The  amplitude  of 
this  undesired  signal  is  reduced  by  using  a  directional 
coupler;  nevertheless  it  represents  a  large  noise  to  the  real 
scattering  signal.  To  resolve  this  problem,  a  standard  noise 
cancellation  technique  is  used.  The  microwave  components 
used  for  noise  cancellation  are  shown  in  the  diagram  (Fig. 
6) .  An  HP  spectrum  analyzer  ( 8569B )  is  used  for  recording 
the  scattering  signal.  It  is  noted  that  the  attenuation  of  the 
directional  coupler  is  frequency  dependent.  Only  test  waves 
with  frequencies  leading  to  more  than  15  dB  attenuation  of 
the  directional  coupler  are  used  in  the  experiment.  Conse- 
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FIG.  5.  Probe  measurement  of  the  plasma  density  distribution  along  the 
direction  transverse  to  the  plasma  layers.  Measurement  is  from  the  central 
point  x  =  0  of  one  layer  to  the  midpoint  x  —  3.24  of  the  next  layer.  The 
minimum  density  102  is  an  estimated  value. 


quently,  the  perturbations  of  the  noise  signal  as  a  result  of  the 
presence  of  the  plasma,  which  in  principle  is  in  the  same 
intensity  level  as  the  scattering  signal,  is  reduced  by  15  or 
more  decibels,  and  will  not  affect  the  measurement  of  the 
Bragg  scattering  signal. 

Presented  in  Fig.  7  are  the  outputs  of  the  spectrum  ana¬ 
lyzers  for  two  cases.  Figure  7(a)  shows  that  no  signal  is 
received  when  there  is  no  plasma.  However,  an  appreciable 
scattering  signal  is  detected,  as  shown  in  Fig.  7(b),  when¬ 
ever  plasma  layers  are  produced.  The  frequency  of  the  test 
wave  is  4.01  GHz,  which  is  much  higher  than  the  cutoff 
frequency.  A  clear  signature  of  Bragg  scattering  has  been 
demonstrated.  The  temporal  evolution  of  the  scattering  sig¬ 
nal  has  also  been  measured.  The  result  for  a  test  wave  with 
frequency  5.5  GHz  is  presented  in  Fig.  8(a).  For  compari¬ 
son,  one  of  the  two  microwave  pulses  (both  are  1.1  /xsec) 
used  for  plasma  generation  is  shown  in  Fig.  8  (b) .  As  we  can 
see,  the  scattering  signal  continues  to  persist  for  about  100 
fisec  after  the  breakdown  pulses  are  turned  off.  This  result 
indicates  that  the  coherent  scattering  process  can  be  very 
effective  even  when  the  plasma  frequency  is  well  below  ( ~  2 
order  of  magnitude)  the  cutoff  frequency  of  the  test  wave. 

The  next  question  is  how  the  experimental  results  com¬ 
pare  with  the  theory.  Considering  a  set  of  AT  parallel  plasma 
slabs  with  thickness  S  and  separation  d,  and  using  the  Bragg 
condition  2d  sin  0  —  nA, ,  where  9  is  the  angle  of  Bragg  scat¬ 
tering,  As  is  the  wavelength  of  the  scattered  wave,  and  n  is  a 
positive  integer,  the  scattering  (reflection)  coefficient  S  of 
the  wave  intensity  is  derived  analytically  and  expressed  as 

S=  |<Fr/e,|2  —  (k2a/2)2[sm(mr5/d)/(nn/d)2]2 

X  [sin  N9,/sm  6,  ]2,  (1) 

where  ef  and  er  are  the  field  amplitudes  of  the  test  wave  and 
scattering  wave,  respectively;  a  —o^/eo,  (co2  +  vj!)1/2;  k, 
and  co,  are  the  wavenumber  and  frequency  of  the  test  wave; 
and  6,  =  (<y,  —  cos)d  sin  9/c  and  cos  =  2irc/As. 

This  reflection  coefficient  is  then  plotted  as  a  function  of 
the  test  wave  frequency,  in  Fig.  9.  By  sweeping  the  test  wave 
frequency,  such  a  dependence  is  also  determined  experimen¬ 
tally  in  a  relatively  small  frequency  range  (4.3— 7.8  GHz) 
and  presented  in  Fig.  9  for  comparison.  The  frequency  de¬ 
pendencies  of  output  intensity  of  the  sweep  generator  and 
the  antenna  gain  of  the  receiving  horn  (S-band  horn  #2) 
have  been  examined  and  taken  into  account  in  calibrating 
the  intensity  of  the  scattering  signals.  Though  a  maximum  of 
eight  layers  can  be  produced,  only  three  of  them  have  signifi¬ 
cant  overlap  along  a  line  of  sight.  Therefore  only  these  layers 
can  significantly  contribute  to  the  Bragg  scattering  process. 
Besides  an  uncalibrated  absolute  magnitude,  the  two  func¬ 
tional  dependencies  are  shown  to  agree  with  each  other  very 
well.  It  is  noted  that  the  separation  d  between  the  two  adja¬ 
cent  plasma  layers  is  related  to  the  wavelength  A0  of  the 
microwave  pulses  and  the  angle  <f>  between  the  propagation 
directions  of  the  two  intersecting  pulses  with  the  relationship 
d  —  A^/ 2  sin  (0/2).  Using  the  Bragg  condition  2d  sin  9 
=  nA,,  the  optimum  frequency  for  Bragg  scattering  is  given 
by 

/«  =  nfo  sin  ( i£/2) /sin  9. 
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FIG.  6.  Microwave  Bragg  scattering  ex¬ 
perimental  setup. 


In  the  present  experiment,  <f>  =  90°  and  6  —  45°,  and 
thus  fs  =  nf0.  This  indicates  that  the  frequency  of  the  test 
wave,  which  satisfies  the  Bragg  condition  for  the  current 
experimental  arrangement,  is  equal  to  the  frequency  and  its 
harmonics  of  the  breakdown  pulses.  Consequently,  the 
breakdown  wave  cannot  be  filtered  out  and  represents  a  very 
strong  noise,  which  prevents  any  meaningful  test  of  Bragg 


(a) 


;tj  , 


(b) 

FIG.  7.  Spectrum  analyzer  CRT  display.  ( a )  No  signal  is  received  when  the 
plasma  is  off.  (b)  The  spectrum  of  scattering  signals  when  plasma  layers  are 
present. 


scattering  at  these  frequencies,  and  in  fact,  also  in  the  neigh¬ 
borhood  frequency  regions.  Although  the  optimum  frequen¬ 
cy  region  for  Bragg  scattering  is  not  examined,  nevertheless 
the  consistency  between  prediction  and  experimental  results 
may  lead  us  to  conclude,  based  on  the  maximum  theoretical 
reflection  coefficient,  that  plasma  layers  can  indeed  be  an 
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(b) 

FIG.  8.  The  time  domain  measurement  of  the  scattering  signal  and  the  mi¬ 
crowave  pulse  used  for  plasma  generation,  10 /isec/division.  (a)  Growth 
and  decay  of  the  scattering  signal  over  a  100/zsec  time  interval,  (b)  A  1.1 
/«ec  pulse  similar  to  that  of  Fig.  2(a). 
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FIG.  9.  The  dependence  of  the  reflectivity  S  of  plasma  layers  on  wave  fre¬ 
quency.  Experimental  and  theoretical  results.  The  0  dB  on  the  vertical  axis 
is  an  arbitrary  reference. 


effective  Bragg  reflector,  especially  if  more  layers  can  be  pro¬ 
duced  for  scattering  purposes. 

IV.  SUMMARY  AND  DISCUSSION 

Plasma  layers  generated  by  two  intersecting  microwave 
pulses  are  used  for  the  study  of  Bragg  scattering.  The  experi¬ 
ment  is  conducted  in  a  large  Plexiglas  chamber  shielded  with 
a  microwave  absorber  so  that  the  microwave  reflection  from 
the  wall  can  be  minimized.  Hence  the  experiment  can  be 
considered  to  be  a  laboratory  simulation  of  conceptualized 
plasma  layers  generated  by  high-power  radio  waves  in  the 
upper  atmosphere,  as  investigated  theoretically  by  Gure¬ 
vich. 

We  first  determine  the  characteristics  of  air  breakdown 
by  powerful  microwave  pulses.  We  have  observed  that  it  is 
very  difficult  to  generate  plasma  in  the  center  of  the  chamber 
with  a  single  pulse.  This  is  mainly  because  the  plasma  gener¬ 
ated  near  the  wall  adjacent  to  the  microwave  horn  causes 
erosion  of  the  tail  of  the  incident  pulse  and  the  pulse  becomes 
too  short,  by  the  time  it  reaches  the  central  region  of  the 
chamber,  to  cause  appreciable  ionization.  However,  this 
problem  is  easily  overcome  when  the  scheme  of  two  inter¬ 
secting  pulses  is  used  for  plasma  generation.  In  this  ap¬ 
proach,  each  pulse  has  its  field  amplitude  below  the  break¬ 
down  threshold  to  avoid  the  tail  erosion.  However,  the  fields 
in  the  intersecting  region  can  add  up  and  exceed  the  break¬ 
down  threshold.  This  scheme  is  most  effective  when  the  two 
pulses  have  the  same  polarization  and  are  coherent.  In  this 
case,  the  wave  fields  form  a  standing  wave  pattern  in  the 
intersecting  region  in  the  direction  perpendicular  to  the  bi¬ 
secting  line  of  the  angle  <f>  between  the  intersecting  pulses. 
Thus  parallel  plasma  layers  with  a  separation 
d  =s  A.f/2  sin  ($/2)  can  be  generated.  This  result  is  shown  in 
Fig.  1. 

Since  there  are  no  electrodes  involved  iv  the  current  ex¬ 
periment  of  air  breakdown,  we  can  determine  the  break¬ 
down  threshold  field  as  a  function  of  the  air  pressure  within 
the  accuracy  of  microwave  probe  measurement.  Two  Pas- 
chen  breakdown  curves  for  the  cases  of  1.1  and  3.3  //sec 


pulses  are  determined  as  shown  in  Fig.  3.  The  appearance  of 
a  Paschen  minimum  can  be  e  xplained  as  the  result  of  break¬ 
down  by  a  short  pulse  that  is  equivalent  to  a  dc  discharge 
with  short  separation  between  electrodes  ( i.e.,  short  electron 
transit  time).  The  result  that  the  breakdown  threshold  by  a 
longer  pulse  (3.3  fisec)  is  lower  than  that  by  a  shorter  pulse 
(1.1  n sec)  agrees  with  the  explanation.  The  characteristic  of 
the  curves  is  also  confirmed  phenomenologically  by  the  var¬ 
ious  degrees  of  tail  erosion  of  the  same  pulse  passing  through 
the  chamber  at  different  pressures,  as  shown  in  Fig.  4. 

An  optical  probe  has  been  used  to  monitor  the  growth 
and  decay  of  airglow  enhanced  by  microwave-generated 
electrons.  Two  processes  are,  in  general,  responsible  for  the 
airglow;  one  is  electron-ion  recombination  and  the  other  is 
impact  excitation  of  natural  gas.  Since  only  weakly  ionized 
plasma  is  generated,  the  second  process  is  believed  to  be 
dominant.  However,  the  second  process  requires  that  the 
electron  energy  exceed  2  eV.  Therefore  the  decay  rate  of 
airglow  intensity  shown  in  Fig.  2(b)  accounts  for  not  only 
the  decay  of  the  electron  density  caused  by  the  dissociative 
attachment  loss  but  also  for  the  decay  of  the  electron  tem¬ 
perature  caused  by  energy  loss  to  the  neutrals  and  the  loss  of 
fast  electrons. 

We  then  conducted  the  Bragg  scattering  experiment 
with  the  produced  plasma  layers.  Both  temporal  evolution  of 
the  scattering  signal  from  a  test  wave  and  the  spectral  de¬ 
pendence  of  the  scattering  coefficient  have  been  examined. 
Good  agreement  between  theoretical  and  experimental  re¬ 
sults  on  the  spectral  dependence  of  the  scattering  coefficient 
has  been  achieved  (Fig.  9).  The  theoretical  result  ( 1 )  shows 
that  the  scattering  coefficient  is  proportional  to  the  square  of 
the  electron  density  and  insensitive  to  the  electron  tempera¬ 
ture.  Therefore  the  electron  decay  rate  after  the  breakdown 
pulses  h^ve  passed  through  can  be  evaluated  from  Fig.  8(a). 
The  figure  shows  that  the  initial  decay  rate  is  about  6  X  104 
sec  “ ' ,  which  is  consistent  with  the  dissociative  attachment 
rate.  It  also  shows  that  the  electron  density  after  70  /rsec  is 
reduced  by  a  factor  of  about  30.  Similarly,  the  electron  decay 
rate  is  also  reduced  over  a  factor  of  10  from  the  initial  decay 
rate.  In  this  region  the  decay  rate  is  consistent  with  the  re¬ 
combination  loss  rate.  This  is  realized  by  the  fact  that  when 
enough  negative  molecule  ions  are  produced  through  the 
electron  attachment  process,  the  detachment  rate  for  elec¬ 
tron  regeneration  is  increased  and  eventually  balances  out 
the  electron  attachment  loss  rate.  Thus  the  dominant  elec¬ 
tron  loss  mechanism  is  shifted  to  the  recombination  process. 
It  should  be  noted  that  the  attachment  process  causes  only 
the  decay  of  free  electrons.  The  ion  density  decays  at  a  rela¬ 
tively  slow  rate  determined  by  the  recombination  process 
and  the  ambipolar  diffusion  process. 
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A  chamber  experiment  is  conducted  to  study  the  propagation  of  high-power  microwave  pulses 
through  the  air.  Two  mechanisms  responsible  for  two  different  degrees  of  tail  erosion  have 
been  identified  experimentally.  The  optimum  pulse  amplitude  for  maximum  energy  transfer 
through  the  air  has  also  been  determined. 


1.  INTRODUCTION 

The  propagation  of  high-power  microwave  pulses 
through  the  atmosphere  has  been  a  subject  with  considerable 
scientific  interest. 1-5  This  is  because  air  breakdown  produces 
ionization  phenomena  that  can  radically  modify  wave  prop¬ 
agation.  Ionization  gives  rise  to  a  space-time-dependent 
plasma  which  attenuates  the  tail  of  the  pulse  but  hardly  af¬ 
fects  the  leading  edge  because  of  the  finite  time  for  the  plas¬ 
ma  to  build  up.  A  mechanism  which  is  called  “tail  erosion” 
plays  the  primary  role  in  limiting  transmission  of  the 
pulse.1-3  Moreover,  the  nonlinear  and  nonlocal  effects 
brought  about  by  the  space-time-dependent  plasma  also  play 
important  roles  in  determining  the  propagation  characteris¬ 
tic  of  the  pulses.3  Therefore,  any  meaningful  theoretical  ef¬ 
fort  requires  a  self-consistent  description  of  the  propagation 
process.  Consequently,  an  experimental  effort  could  be  more 
relevant  and  useful. 

Basically,  there  are  two  fundamental  issues  to  be  ad¬ 
dressed.  One  concerns  the  optimum  pulse  characteristics  for 
maximum  energy  transfer  through  the  atmosphere  by  the 
pulse.  The  second  concern  is  maximizing  the  ionizations  in 
the  plasma  trail  following  the  pulse.  In  general,  these  two 
concerns  are  interrelated  and  must  be  considered  together. 
This  is  because  in  order  to  minimize  the  energy  loss  in  the 
pulse  before  reaching  the  destination,  one  has  to  prevent  the 
occurrence  of  excessive  ionization  in  the  background  air. 
Otherwise,  the  overdense  plasma  can  cutoff  the  propagation 
of  the  remaining  part  of  the  pulse  and  cause  the  tail  of  the 
pulse  to  be  eroded  via  the  reflection  process.  This  process  is 
believed  to  be  far  more  severe  in  causing  tail  erosion  than  the 
normal  process  attributed  to  ionization  and  heating.  Once 
this  process  occurs,  the  remaining  pulse  will  become  too  nar¬ 
row  to  ionize  dense  enough  plasma  whose  nonlinear  effect  is 
thought  to  help  in  sufficient  focusing  to  compensate  for  wave 
spreading  beyond  the  Fresnel  distance. 

The  purpose  of  our  experimental  effort  is  to  understand 
the  fundamental  behavior  of  tail  erosion  and  address  the 
question  of  how  energy  loss  depends  on  basic  parameters 
such  as  pulse  intensity  and  width,  and  background  pressure. 
The  experimental  data  can  then  be  incorporated  for  the  de¬ 
velopment  of  a  useful  theoretical  model  for  a  self-consistent 
derivation  of  pulse  propagation. 


II.  TAIL  EROSION  MECHANISMS 

Experiments  have  been  conducted  in  a  large  chamber 
made  of  a  2-ft  cube  of  plexiglass  and  filled  with  dry  air  at 
various  pressures.  A  microwave  pulse  is  fed  into  the  cube  by 
an  5-band  microwave  horn  placed  at  one  side  of  the 
chamber.  A  second  5-band  horn  placed  at  the  opposite  side 
of  the  chamber  is  used  to  receive  the  transmitted  pulse.  The 
chamber  is  shielded  with  microwave  absorbers  so  that  the 
microwave  reflection  from  the  walls  and  nearby  structures 
can  be  minimized.  The  microwave  power  is  generated  by  a 
single  magnetron  tube  (OKH  1448)  driven  by  a  soft-tube 
modulator.  The  magnetron  produces  1  MW  peak  output 
power  at  a  frequency  of  3.27  GHz.  The  modulator  uses  a 
pulse-forming  network  (PFN)  having  a  pulse  width  of  1.1 
fj.s  and  a  repetition  rate  of  60  Hz.  The  PFN  can  also  be  ad¬ 
justed  to  increase  the  pulse  width  to  3.3  /xs,  while  reducing 
the  repetition  rate  to  20  Hz.  Figure  1  shows  the  block  dia¬ 
gram  of  the  experimental  setup. 

Breakdown  of  air  was  detected  either  visually,  as  the 
first  sign  of  a  glow  in  the  chamber,  or  as  the  distortion  in  the 
shape  of  the  pulse  received  by  the  horn  placed  at  the  opposite 
side  of  the  chamber.  The  dependence  of  the  breakdown 
threshold  field  on  the  pressure  can  thus  be  manifested  by  a 
dependence  of  the  degree  of  distortion  of  the  tail  portion  of 
the  transmitted  pulse  through  the  chamber.  A  series  of  snap 
shots  demonstrating  this  behavior  is  presented  in  Fig.  2.  In 


screen  room 


FIG.  1.  Experimental  setup. 
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FIG.  2.  Pulses  propagating  through  chamber  at  different  air  pressures,  (a)  p  =  8,  (b)  p  =  4,  (c)  p  -  1,  (d)  p  -  0.2,  (e)  p  =  0.07,  and  (f)  p  =  0.05  Torr. 
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the  high-pressure  regions  pictured,  the  breakdown  threshold 
is  expected  to  be  high  and,  therefore,  very  little  ionization 
can  occur  and  the  pulse  can  pass  through  the  chamber  al¬ 
most  without  any  loss  (or  distortion).  However,  as  the  pres¬ 
sure  drops,  the  breakdown  threshold  also  decreases  before 
reaching  the  minimum,  hence,  more  ionization  occurs  and 
also  more  distortion  to  the  pulse.  The  distortion  always 
starts  from  the  tail  of  the  pulse  (i.e.,  tail  erosion)  because  it 
takes  finite  time  for  the  plasma  to  build  up  and,  thus,  maxi¬ 
mum  absorption  of  pulse  energy  by  the  generated  electrons 
always  appears  in  the  tail  of  the  pulse.  Consequently,  the 
leading  edge  of  the  pulse  is  usually  not  affected.  Between  1 
and  2  Torr,  the  pulse  appears  to  suffer  maximum  tail  erosion 
and  hence  only  very  narrow  leading  edge  can  pass  through 
the  chamber.  The  tail  erosion  becomes  weak  again  for  a 
further  decrease  in  the  pressure  and  can  eventually  be  pre¬ 
vented  once  the  pressure  becomes  so  low  that  the  breakdown 


threshold  power  exceeds  the  peak  power  of  the  incident 
pulse. 

Tail  erosion  is  a  common  phenomenon  appearing  in  the 
propagation  of  a  high-power  microwave  pulse  (HMP).  This 
phenomenon  is  further  demonstrated  by  the  snap  shots  pre¬ 
sented  in  Fig.  3,  where  1.1  /zs  pulses,  with  four  consecutively 
increasing  amplitudes,  are  transmitted  into  the  chamber  of  1 
Torr  pressure  from  one  side  and  received  at  the  opposite 
side.  The  first  pulse  has  amplitude  below  the  breakdown 
threshold,  and  hence,  nothing  is  expected  to  happen.  Conse¬ 
quently,  the  received  pulse  shape,  as  shown  in  Fig.  3(a),  is 
undistorted  from  that  of  transmitted  pulse.  Once  the  ampli¬ 
tude  exceeds  the  breakdown  threshold,  more  tail  erosion  oc¬ 
curred  to  the  larger  amplitude  pulses,  as  is  observed  by  the 
subsequent  three  snapshots  [Figs.  3(b)— 3(d) ].  This  is  be¬ 
cause  the  increase  of  the  ionization  rate  with  field  amplitude 
allows  more  electrons,  which  attenuate  the  pulse,  to  build 
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up.  Now  let  us  focus  on  the  last  two  picutures  [Figs.  3(c)  caused  through  reflection  by  the  self-generated  overdense 

and  3(d)  ] .  Pulses  have  been  eroded  strongly  in  both  cases.  plasma  screen.  These  two  processes  are  also  verified  by  the 

However,  a  clear  distinction  between  the  two  cases  is  no-  reflected  power  level  measured  for  each  case.  As  shown  in 

ticed.  In  one  case  corresponding  to  the  third  picture  [Fig.  Fig.  4,  the  snap  shots  of  Figs.  4(b)  and  4(d)  are  the  reflected 

3(c)  ],  the  erosion  to  the  tail  of  the  pulse  is  not  complete.  In  pulse  shape  corresponding  to  each  received  pulse  of  Figs, 

other  words,  the  received  pulse  width  extends  to  the  original  4(a)  and  4(c),  respectively.  As  shown  by  the  last  set  of  pic- 

width.  In  the  other  case  [Fig.  3(d)  ],  a  large  portion  of  the  tures  [Figs.  4(c)  and  4(d)  ],  strong  reflection  and  complete 

pulse  is  more  or  less  eroded  completely  during  the  finite  erosion  are  observed  consistently. 

propagation  period.  Obviously  it  is  a  different  mechanism  In  order  to  avoid  cutoff  reflection,  the  power  of  the  pulse 

responsible  for  the  second  case.  The  ionization  frequency  should  be  lower  than  a  critical  power  Pr  which  is  defined  as 

becomes  so  large  in  the  second  case  that  the  electron  density  the  minimum  required  power  for  generating  an  overdense 

exceeds  the  cutoff  density  of  the  wave  before  the  whole  pulse  plasma  screen.  When  an  overdense  plasma  screen  is  formed 

passes  through.  The  overdense  plasma  screen  reflects  the  the  shape  of  the  reflected  pulse  changes  drastically  and  can 

remaining  portion  of  the  pulse  and  causes  even  more  severe  be  monitored  easily.  Thus,  the  critical  power  Pr  can  be  deter- 

tail  erosion.  In  summary,  two  mechanisms  responsible  for  mined.  This  critical  power  varies  in  general  with  the  pres- 

the  tail  erosion  are  identified.  One  is  due  to  attenuation  by  sure.  Measurements  are  made  to  determine  this  functional 

the  self-generated  underdense  plasma.  The  other  one  is  dependency  as  shown  in  Fig.  5,  where  P  =  Pr/Pc  is  the  criti- 


(b)  (d) 

FIG.  4.  Received  pulses  and  the  corresponding  reflected  pulses,  (a)  Envelope  of  received  pulse,  (b)  Reflected  pulse  corresponding  tocase(a).  (c)  Received 
pulse  with  larger  incident  amplitude  than  that  of  case  (a),  (d)  Reflected  pulse  corresponding  to  case  (c). 
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FIG.  5.  The  dependence  of  the  normalized  critical  power  P  =  Pr/Pc  on 
pressure  for  1. 1  /rs  pulse. 


cal  power  normalized  to  the  breakdown  threshold  power  Pc 
which  has  been  presented  elsewhere.'’  Due  to  the  limited 
available  microwave  power,  only  the  region  0.2-10  Torr  is 
examined.  Outside  this  pressure  region,  the  breakdown 
threshold  power  becomes  too  high  to  be  exceeded  by  the 
maximum  output  power  of  our  microwave  source.  Never¬ 
theless,  this  is  the  pressure  region  of  main  concern  because 
air  breakdown  has  its  lowest  threshold. 

III.  DISCUSSION 

Our  study  indicates  that  an  increase  of  pulse  amplitude 
may  not  help  to  increase  the  energy  transfer  by  the  pulse. 
This  is  because  two  tail  erosion  mechanisms  are  at  play  to 
degrade  the  energy  transfer.  A  demonstration  is  presented  in 
Fig.  6,  where  the  growth  and  decay  of  airglow  enhanced  by 
electrons  through  air  breakdown  by  3.3  /us  pulse  are  record¬ 
ed  for  two  different  power  levels.  In  Fig.  6  ( a )  the  power  level 
is  below  the  critical  power  and  the  airglow  grows  for  the 
entire  3.3  /us  period  of  the  initial  pulse  width.  As  power  is 
increased  beyond  the  critical  value,  the  initial  growth  of  the 
airglow  becomes  faster  as  shown  in  Fig.  6(b).  However,  it  is 
also  shown  in  Fig.  6(b)  that  the  airglow  saturates  at  about 
the  same  level  as  that  of  Fig.  6(a).  Moreover,  the  airglow 
already  starts  to  decay  even  before  the  3.3  /us  period.  In  other 
words,  cutoff  reflection  happening  in  the  second  case  limits 
the  energy  transfer  by  the  pulse.  The  additional  energy  add¬ 
ed  to  the  pulse  is  wasted  by  reflection.  The  way  to  solve  the 
problem  is  either  to  lower  the  amplitude  of  the  pulse  or  to 
narrow  the  pulse  width  so  that  the  propagation  loss  can  be 
minimized.  The  remaining  question  is  about  the  optimum 
parameters  for  the  pulse  in  order  to  achieve  effective  ioniza¬ 
tion  for  a  self-focusing  purpose.  Unfortunately,  this  question 
cannot  easily  be  answered  by  the  chamber  experiments.  This 
is  because  the  effect  of  pressure  gradient  and  large  propaga¬ 
tion  distance  cannot  be  incorporated  in  the  investigation. 
Nevertheless,  the  understanding  of  the  fundamental  behav- 


(b) 


FIG.  6.  The  growth  decay  of  enhanced  airglow  (upper  trace).  Pulses  of  3.3 
fis  duration  and  at  two  different  power  levels  P,  (for  (a)  |  and  P,  (for  (b)  ] 
are  used  for  causing  air  breakdown,  where  P-JP ,  =  1 .32.  The  lower  trace  of 
each  photo  represents  the  time  dependence  of  the  amplitude  of  the  micro- 
wave  pulse.  The  horizontal  scales  are  1  /js/div. 

ior  of  tail  erosion  and  the  set  of  experimental  data  obtained 
should  be  very  useful  for  guiding  the  development  of  a  prac¬ 
tical  theoretical  model. 
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Two  crossed  microwave  pulses  are  used  to  generate  plasma  in  their  intersection  region  inside  a 
chamber.  It  is  shown  that  the  frequency  of  each  output  pulse  is  up-shifted.  The  percentage  of  frequency 
shift  is  consistent  with  the  theory  of  Wilks,  Dawson,  and  Mori.  Besides  frequency  shift  and  some  at¬ 
tenuation,  the  similarity  between  the  spectral  distributions  of  the  output  and  input  pulses  indicates  that 
there  is  very  little  distortion  to  the  pulse  shape. 

PACS  numbers:  52.40.Db 


Recently,  Wilks,  Dawson,  and  Mori1  have  studied  an 
initial-value  problem  concerning  the  propagation  of  an 
electromagnetic  wave  through  a  gas  which  is  ionized  in  a 
very  short  time  interval  on  the  order  of  the  period  of  this 
wave.  The  results  of  their  ID  and  2D  particle-in-cell 
computer  simulations  show  that  the  frequency  of  the 
wave  is  up  shifted  and  confirm  their  linear  theory  for  the 
ideal  case  of  instantaneous  plasma  creation.  The  physi¬ 
cal  mechanism  of  this  frequency  up-conversion  phenom¬ 
enon  becomes  quite  clear  by  noting  that  the  wavelength 
of  the  wave  does  not  change  after  the  ionization  of  the 
gas  since  the  wave  does  not  encounter  any  spatial  varia¬ 
tion  in  the  medium;  however,  the  sudden  reduction  of  the 
index  of  refraction  of  the  gas  by  converting  it  into  a  plas¬ 
ma  forces  the  wave  to  propagate  with  larger  phase  veloc¬ 
ity  and,  thus,  to  oscillate  with  higher  frequency  at  subse¬ 
quent  times. 

However,  using  a  single  high-power  laser  pulse  for 
plasma  generation  and,  subsequently,  for  its  frequency 
up-conversion  purpose  was  shown  to  be  unpractical  in 
the  previous  experiment  by  Yablonovitch.2  His  experi¬ 
ment  showed  only  spectral  broadening  in  the  CO2  laser 
pulse  transmitting  through  a  rapidly  growing  plasma. 
We  believe  that  this  is  because  plasma  is  preferentially 
generated  at  the  boundary  where  input  pulse  and  ionized 
medium  first  encounter. 

In  this  work,  we  use  two  crossed  microwave  beams 
with  parallel  polarization  for  establishing  an  initial-value 
problem.  In  doing  so,  we  can  adjust  the  intensity  of  each 
beam  so  that  a  discharge  of  background  gas  occurs  only 
in  the  intersection  region  of  the  two  beams.  The  experi¬ 
ment  is  conducted  in  a  chamber  of  a  2-ft  Plexiglas  cube. 
The  carrier  frequency  of  the  1.1 -//sec  microwave  pulses 
used  in  the  experiment  is  3.27  GHz.  Shown  in  Fig.  1(a) 
is  the  envelope  of  each  pulse.  The  chamber  is  hence 
more  than  six  wavelengths  long  and  is  long  enough  to 
shield  the  boundary  effect  from  the  wave  over  a  single 
wave  period.  The  microwave  power  is  generated  by  a 
single  magnetron  tube  (OKH1448)  driven  by  a  soft  tube 
modulator  having  a  repetition  rate  of  60  Hz.  The  peak 
output  power  of  the  tube  is  1  MW.  Since  we  are  in¬ 


terested  in  an  initial-value  problem,  we  have  to  ionize  the 
gas  simultaneously  over  the  chamber.  However,  using  a 
single  pulse  can  only  generate  a  localized  plasma  near 
the  chamber  wall.  A  second  pulse  provided  by  the  same 
magnetron  tube  is  fed  into  the  cube  through  a  second  S- 
band  microwave  horn  placed  at  a  right  angle  to  the  first 
one.  With  such  an  arrangement,  the  power  of  each  pulse 
can  be  reduced  to  below  the  breakdown  threshold  of  the 
gas  and  breakdown  of  gas  can  occur  only  in  the  central 


<b) 
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FIG.  I.  (a)  Envelope  of  a  I  -//sec  microwave  pulse,  (b)  A 
photo  of  plasma  layers  generated  by  two  cross  pulses  having 
parallel  polarization. 
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region  of  the  chamber  where  the  two  pulses  intersect. 
The  undesired  tail  erosion  to  the  pulses  by  the  generated 
plasma  can  also  be  minimized.  The  wave  fields  add  up 
to  form  a  standing-wave  pattern  in  the  intersecting  re¬ 
gion  in  the  direction  perpendicular  to  the  bisecting  line 
of  the  angle  between  the  two  intersecting  pulses.  Thus, 
parallel  plasma  layers  with  a  separation  d—Xjyfi can  be 
generated,  where  Xo  is  the  wavelength  of  the  wave.  Two 
phase  shifters  are  used  to  extend  the  region  covered  by 
the  plasma  layers.  A  maximum  of  eight  layers  detected 
visually  can  be  generated  to  cover  almost  the  entire  path 
of  the  first  pulse  in  the  chamber. 3  A  photo  of  the  plasma 
layers  is  presented  in  Fig.  1  (b),  where  only  seven  layers 
are  clearly  shown.  A  block  diagram  of  the  experimental 
setup  is  shown  in  Fig.  2. 

Although  the  produced  plasma  is  not  spatially  uni¬ 
form,  it  extends  over  several  wavelengths  and  has  a  tem¬ 
poral  growing  average  density  along  the  path  of  the  first 
pulse.  The  periodic  spatial  variation  of  plasma  density 
may  cause  scattering  and,  hence,  additional  attenuation 
to  the  incident  pulse.  Nevertheless,  it  should  not  affect 
the  observation  on  frequency  up-shift  phenomenon. 
Since  our  pulse  length  is  about  300  m  long,  which  is 
much  longer  than  the  length  0.6  m  of  the  chamber,  one 
may  divide  the  pulse  into  500  segments  and  each  seg¬ 
ment  of  the  pulse  experiences  a  frequency  shift  while  it  is 
inside  the  chamber.  Therefore,  the  maximum  frequency 
shift  on  the  output  pulse  is  just  the  maximum  frequency 
shift  on  each  segment.  Consequently,  only  a  small  frac¬ 
tion  of  total  density  change  over  the  pulse  period  can 
contribute  to  the  frequency  shift.  In  order  to  enhance 
frequency  shift,  one  would  have  to  increase  the  ioniza¬ 
tion  rate.  However,  this  cannot  be  done  by  simply  in¬ 
creasing  the  pulse  amplitude  because  gas  breakdown  can 
occur  locally  near  the  chamber  wall,  and  the  nature  of 
initial-value  problem  is  demolished.  We  find  that  if  we 
mix  dry  air  with  argon  gas  with  total  pressure  about  2.3 
Torr,  maximum  frequency  up-shift  is  achieved.  This  is 
because  dry  air  has  a  high  breakdown  threshold  which 


enables  us  to  increase  the  pulse  amplitude  without  caus¬ 
ing  serious  tail  erosion  and  argon  is  helping  the  initial 
ionization  for  the  subsequent  electron-density  growth. 
Shown  in  Fig.  3  is  the  frequency  spectra  of  output  pulses 
with  and  without  plasma  generation.  The  corresponding 
input  spectra  were  also  examined.  No  noticeable  dif¬ 
ference  between  the  two  input  spectra  was  observed.  A 
1.5-MHz  frequency  up-shift  is  clearly  demonstrated  in 
Fig.  3.  It  is  noted  that  the  sweep  time  of  the  spectrum 
analyzer  is  50  sec.  Thus,  each  recorded  spectrum  of  Fig. 
3  represents  an  average  sampling  spectrum  of  3000 
pulses.  If  the  frequency  shift  of  each  pulse  is  not  the 
same,  the  recorded  spectrum  should  reveal  a  mean  fre¬ 
quency  shift  as  well  as  spectral  broadening.  However,  no 
noticeable  spectral  broadening  appeared  in  the  recorded 
spectrum  (curve  2  of  Fig.  3),  and  no  noticeable  change 
in  the  amount  of  frequency  up-shift  over  many  sweeping 
periods  can  be  detected  either.  Besides  a  frequency  shift 
and  an  attenuation  caused  by  the  plasma  generation,  the 
similarity  between  the  spectrum  of  the  two  pulses  indi¬ 
cates  that  the  distortion  to  the  pulse  shape  due  to  the 
plasma  generation  is  very  small. 

In  the  density  range  of  interest  (<  10 11  cm-3),  the 
dominant  loss  mechanism  of  free  electrons  in  the  dry  air 
is  through  their  attachment  with  neutral  molecules.  The 
maximum  attachment  rate4  va  of  dry  air  at  2.3  Torr  is 
about  3.22xl05  sec-1.  Hence,  breakdown  occurs  only 
when  the  ioniza*ion  rate  u,  is  larger  than  the  attachment 
rate  va.  The  relationship  between  u,  and  va  is  shown4  to 
be  Vi/va~(El Ec )5  3  for  a  long  pulse  and  1.3  <E/EC 
<  3.5,  where  Ec  is  the  breakdown  threshold  field  of  the 
dry  air  and  can  be  expressed4  as 

£cA28(Am/2.7x  10I9)(1  +a>$/v}) 1/2  kV/cm  , 

Nm  is  the  number  of  neutrals  per  cubic  centimeter,  and 


FIG.  3.  Frequency  spectra  of  two  output  pulses  with  same 
input.  Curve  1  is  the  spectrum  of  the  pulse  passing  through 
the  chamber  (at  the  pressure  ~  I  atm)  without  plasma  genera¬ 
tion  and  curve  2  represents  the  spectrum  of  the  pulse  passing 
through  a  rapidly  growing  plasma  generated  inside  the 
chamber  (at  the  pressure  ~2.3  Torr). 
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vc  is  the  electron-neutral  collision  frequency.  In  this  ex¬ 
periment,  (0&/v}^>  1,  thus  Ec—o> o/vc.  Since  vc  is  pro¬ 
portional  to  electron  speed  and  in  the  present  case  the 
electron  quiver  speed  eE/mcoo  is  much  larger  than  the 
electron  thermal  speed,  then  vc~E  and  this  leads  to 
v,-— (£/£c)53~£ l0,6;  a  very  strong  dependence  of  the 
ionization  frequency  v,  on  £  is  shown.  It  is  a  disadvan¬ 
tageous  feature  for  controlling  the  frequency  shift  by 
self-generated  plasmas.  On  the  other  hand,  it  becomes 
an  advantageous  feature  for  possibly  minimizing  the 
spectral  broadening  on  the  frequency-shifted  output 
pulse  as  indicated  by  the  result  shown  in  Fig.  3.  The 
amount  of  electron-density  change  contributed  to  fre¬ 
quency  shift  of  each  pulse  segment  depends  not  only  on 
the  ionization  frequency  v,  but  also  on  the  available  time 
period  t  during  which  the  initial  effect  dominates  over 
the  boundary  effect,  where  r  is  a  time-dependent  func¬ 
tion.  In  the  beginning  of  the  pulse,  the  background  elec¬ 
tron  density  is  low  and  hence  the  boundary  effect  is 
weak.  Thus,  the  initial  effect  can  prevail  over  the  bound¬ 
ary  effect  for  almost  the  entire  period  of  the  pulse  seg¬ 
ment;  i.e.,  ro“//c,  the  transit  time  of  the  pulse  through 
the  chamber,  r  decreases,  however,  as  the  electron  den¬ 
sity  increases.  The  density  change  of  plasma  electrons 
available  for  the  frequency  shift  of  each  pulse  segment  is 
then  given  by  An(/no~v,r.  Since  each  segment  of  the 
pulse  suddenly  experiences  half  of  such  a  fractional 
electron-density  increase  on  the  average,  while  passing 
through  the  chamber,  the  dispersion  relation  of  the  car¬ 
rier  of  the  pulse  should  also  change  accordingly  in  time 
from  arf—ioto+kfic2  to 

to}  —  (tOpo+AtOp )2  +  k&c2“tofi  +  2o)poAo)p + Aa } . 

In  the  region  Atop  <&topo,  we  have 

rn/ss  <u0(l  +2(OpoAa)p/a>$) 1/2 o>o(l  +topo Atop/a>$) . 

Since  AtOp/(Op0  “  t  (Ane/n0)  —  |  An,/no,  we  then  have  the 
frequency  up-shift 

A/0//0-  j  (An,/no)<i>po/W— v/T/io . 

In  general,  this  frequency  shift  is  a  function  of  time.  It 
would  be  the  case  of  a  rectangular  pulse  and  lead  to  a 
spectral  broadening  because  each  pulse  segment  experi¬ 
ences  a  slightly  different  frequency  shift.  However,  no 
apparently  spectral  broadening  is  shown  in  the  result 
(curve  2)  presented  in  Fig.  3.  On  the  other  hand,  the 
pulse  envelope  shown  in  Fig.  1(a)  is  not  exactly  a  rec¬ 
tangular  pulse  either.  It  shows  that  the  pulse  has  a  de¬ 
creasing  amplitude  (on  average).  It  means  that  v,  de¬ 
creases  while  no  increases.  It  is  verified  by  the  Bragg 
scattering  result.  Since  the  generated  plasma  is  a  set  of 
parallel  layers,  one  can  use  the  Bragg  scattering  method 3 
to  monitor  the  growth  and  decay  of  plasma  electrons, 
where  the  scattering  coefficient  is  proportional  to  rtfi(t). 
As  shown  in  Fig.  4,  the  intensity  of  the  scattering  signal 


FIG.  4.  The  growth  and  decay  of  the  intensity  of  Bragg 
scattering  signal  during  the  first  2-psec  period. 


after  emerging  out  of  the  noise  level  increases  with  a  de¬ 
creasing  rate.  In  order  to  minimize  the  spectral  broaden¬ 
ing  on  the  frequency-shifted  output  pulse,  it  requires  that 
v,  r/io~ const  over  a  large  portion  of  pulse  length.  Con¬ 
sidering  the  portion  of  the  pulse  from  its  peak  to  the 
half-peak  point  (— 0.7-/tsec  interval),  a  self-consistent 
set  of  the  solution  of  the  rate  equation  dno/d/*v,n0 
given  by  n0-«i(l  + v,o//1.25)125,  v,  -  vi0/(  1  +  vi0t/ 
1.25),  and  t"to/(1  +  vi0//1.25)025,  with  /ii •“2.5x  107 
cm-3,  v, 0 "■  2.94 x  109  sec _l,  and  t0”*//c~ 2  nsec,  leads 
to  a  constant  frequency  shift  with  a  value  about  1.4 
MHz  consistent  with  the  experimental  result.  Using  the 
function  v,(f)  and  the  relation  v,  —£106,  the  envelope  of 
a  0.7-psec  pulse  can  be  reconstructed  and  is  found  to 
overlap  closely  with  the  original  pulse  envelope  in  the 
corresponding  interval  [Fig.  1(a)  with  average  ampli¬ 
tude].  This  explains  why  there  is  no  apparently  spectral 
broadening  appearing  in  curve  2  of  Fig.  3.  It  is  noted 
that  at  the  end  of  the  pulse,  the  average  plasma  density 
is  about  the  cutoff  density  of  each  input  pulse,  i.e., 

{(Op)  ~  (0$  “  J  (Op  I  peak- 

In  summary,  we  have  demonstrated  experimentally 
that  it  is  possible  to  up-shift  the  carrier  frequency  of  an 
electromagnetic  pulse  by  passing  it  through  a  rapidly 
growing  plasma.  Though  the  achieved  frequency  shift  is 
small,  it  is  mainly  because  the  pulse  length  of  the  present 
experiment  is  much  longer  than  the  size  of  generated 
plasma.  If  the  pulse  length  can  be  reduced  to  be  compa¬ 
rable  to  the  length  of  the  chamber  (e.g.,  using  2-nsec 
pulses),  the  field  intensity  of  each  pulse  can  be  increased 
considerably  without  causing  gas  breakdown  near  the 
chamber  wall.  Thus,  thfc  ionization  rate  of  gas  in  the  in¬ 
tersecting  regions  of  two  pulses  can  be  increased  by  a 
factor  of  about  500,  and  a  frequency  shift  A/o//o~0.26 
is  expected.  In  order  to  minimize  spectral  broadening  on 
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output  pulses,  a  proper  shaping  of  the  input  pulses  is  re¬ 
quired. 
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A  chamber  experiment  is  conducted  to  study  the  propagation  of  a  high-power  microwave  pulse.  The  results  show  that  the  pulse 
is  experiencing  frequency  up-shift  while  ionizing  the  background  air  if  the  initial  carrier  frequency  of  the  pulse  is  higher  than  the 
electron  plasma  frequency  at  the  incident  boundary.  Such  a  frequency  autoconversion  process  may  lead  to  reflectionless  propa¬ 
gation  of  a  high-power  microwave  pulse  through  the  atmosphere. 


1.  Introduction 

During  the  past  decade,  there  has  been  a'  consid¬ 
erable  effort  in  the  study  of  the  propagation  of  high 
power  microwave  pulses  through  the  atmosphere  for 
its  applications  in  high  power  radar  and  directed  en¬ 
ergy  systems  [1-6].  This  study  has  also  been  a  sub¬ 
ject  of  considerable  scientific  interest  because  air 
breakdown  produces  ionization  phenomena  tha,  can 
radically  modify  wave  propagation  [7],  Both  fluid 
and  kinetic  transport  codes  have  been  developed  to 
give  a  self-consistent  description  of  the  propagation 
process.  Many  laboratory  experiments  have  also  been 
conducted  to  verify  the  results  of  computer  codes. 

After  air  breakdown,  a  space-time  dependent 
plasma  is  generated.  It  attenuates  the  tail  of  the  pulse 
but  hardly  affects  its  leading  edge  because  it  takes 
finite  time  for  the  plasma  to  build  up.  A  phenome¬ 
non  which  is  called  “tail  erosion”  plays  the  primary 
role  in  limiting  the  transmission  of  pulse  energy.  Two 
mechanisms  responsible  for  the  tail  erosion  phenom¬ 
enon  have  been  identified  experimentally  in  our  pre¬ 
vious  work  [7],  The  first  one  is  due  to  attenuation 
by  the  self-generated  underdense  plasma.  It  is  shown 
that  the  erosion  to  the  tail  of  the  pulse  is  not  com¬ 
plete  and  the  loss  of  the  pulse’s  energy  may  still  be 
within  the  tolerable  level.  The  second  one  is  caused 
through  reflection  by  the  self-generated  overdense 
plasma  screen.  In  this  case,  as  shown,  a  large  portion 


of  the  pulse  is  more  or  less  eroded  away  during  the 
finite  propagation  period.  The  second  mechanism  is 
far  more  severe  in  causing  energy  loss  of  the  pulse 
than  the  first  one.  Once  the  reflection  mechanism  is 
in  play,  the  remaining  pulse  will  become  too  narrow 
to  carry  enough  energy  as  desired. 

One  easy  way  to  avoid  cutoff  reflection  is  to  re¬ 
strict  the  power  of  the  pulse  with  a  given  pulse  length 
to  be  lower  than  a  critical  power  level.  However,  this 
critical  power  level  is  usually  too  low  for  many  prac¬ 
tical  applications.  Therefore,  it  becomes  essential  to 
explore  means  for  more  powerful  electromagnetic 
pulses  to  achieve  relfectionless  propagation  through 
the  atmosphere.  Recently,  Gildenburg  et  al.  [  8  ]  have 
pointed  out  that  it  is  indeed  possible  to  achieve  re¬ 
flectionless  propagation  of  powerful  EM  pulses  in  an 
ionizing  medium.  Their  analytical  and  numerical  re¬ 
sults  [  9  ]  show  that  a  frequency  autoconversion  pro¬ 
cess  is  in  play  in  a  nonstationary  plasma.  In  partic¬ 
ular,  when  the  ionization  frequency  of  the 
background  air  becomes  much  larger  than  the  elec¬ 
tron-neutral  collision  frequency  v,  the  carrier  fre¬ 
quency  m  of  the  pulse  will  always  shift  upward  with 
the  growth  of  plasma  frequency  so  that  the  plasma 
during  breakdown  will  always  remain  transparent  to 
the  pulse  (i.e.  cu>wpe). 

In  this  paper,  we  describe  a  chamber  experiment 
conducted  to  demonstrate  the  frequency  autocon¬ 
version  process.  The  chamber  is  made  of  a  2  foot  cube 
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of  plexiglass.  A  microwave  pulse  is  fed  into  the  cube 
by  an  S-band  microwave  horn  placed  at  one  side  of 
the  chamber.  A  second  S-band  horn  placed  at  the  op¬ 
posite  side  of  the  chamber  is  used  to  receive  the 
transmitted  pulse.  The  spectra  of  the  incident  pulse 
and  transmitted  pulse  are  then  compared.  The  mi¬ 
crowave  power  is  generated  by  a  single  magnetron 
tube  (OKH1448)  driven  by  a  soft  tube  modulator. 
The  magnetron  produces  1  MW  peak  output  power 
at  a  frequency  of  3.27  GHz.  The  modulator  uses  a 
pulse  forming  network  (PFN)  having  a  pulse  width 
of  1  ps  and  a  repetition  rate  of  60  Hz. 

In  section  2,  a  theoretical  model  describing  pulse 
propagation  in  an  ionizing,  air  is  introduced.  The 
equation  describing  the  frequency  autoconversion 
process  is  also  derived  and  discussed.  The  pulse 
propagation  experiment  is  presented  in  section  3  in 
which  the  experimental  results  are  also  included.  A 
discussion  of  the  work  is  given  in  section  4. 


2.  Modeling  of  the  propagation  of  a  powerful  EM 
pulse  in  ionizing  air 


The  processes  which  are  relevant  to  the  plasma 
formation  include  ionization,  attachment,  recombi¬ 
nation,  and  diffusion.  When  the  power  level  of  the 
pulse  in  consideration  far  exceeds  the  breakdown 
threshold  of  the  air,  however,  the  ionization  rate  v{ 
of  the  air  dominates  over  all  the  loss  rates  of  the  in¬ 
duced  plasma.  Thus,  the  rate  equation  for  the  elec¬ 
tron  density  during  the  pulse  period  can  be  simpli¬ 
fied  by 


where  w  =4jtnceJ/me  is  linearly  proportional  to  the 
electron  density  n,  and  v,  is  proportional  to  the  in¬ 
tensity  of  the  pulse. 

The  propagation  of  an  EM  pulse  in  such  an  ion¬ 
izing  medium  is  governed  by  the  wave  equation 


where  the  electron  momentum  equation  has  been 
used  to  determine  the  current  density  induced  by  the 


wave  field  in  the  plasma;  v  is  the  electron-neutral 
collision  frequency  and  v=  v+  v,. 

Eq.  (2)  is  analyzed  by  letting  1 

£(jc,  t)=A(x,  t) 

and  assuming  that 

I1"-4  *  i* 


and 


for  a  pulse  which  contains  many  oscillations.  Fur¬ 
thermore,  the  local  frequency  aj=d<p/dt  and  wave- 
number  k=—d<t>/dx  are  defined.  Thus,  (2)  can  be 
reduced  approximately  to 


^-cjA2  +  ~~  kc2A2  =  — 
dt  dx 


ioI'VcdA2 

w2+v2 


(3) 


and 

”2-k2c2=<{'-Jh'2)-  (4) 

From  the  definition  of  oj  and  k,  it  is  required  that 


dx  dl 


(5) 


Taking  the  partial  time  derivative  of  (4)  and  us¬ 
ing  relation  (5),  the  equation  describing  the  fre¬ 
quency  autoconversion  is  derived  to  be 


d  ,  d  1  -vv/(a)2  +  v2)  d  , 

co  —  ,  .,  -it  ■>  ,  -2 ,  ■*  ,  tu  pc 

dt  1  dx  1  -a)itPv/(ar+v-)~  dt 

(6) 


where  r,=3c o/dk  is  the  local  group  velocity  of  the 
pulse. 

Under  the  conditions  that  vv<&q)2  +  v2  and 
wje  <cj2+v2,  (6)  reduces  to  the  equation  derived 
by  Gildenburg  et  al.  [8], 


(7) 


where  d/d i=d/dt+vtd/dx  is  the  total  time  deriva¬ 
tive  along  the  group  trajectory  of  the  pulse. 

The  implication  of  cq.  ( 7 )  is  illustrated  as  follows. 
Consider  a  rectangular  pulse  as  shown  in  fig.  1.  the 
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Fig.  I .  Envelope  of  a  rectangular  pulse  and  the  pulse-generated 
electron  density  distribution. 

density  of  the  induced  plasma  at  a  given  time  t  is  ex¬ 
pected  to  grow  spatially  from  the  background  level 
at  the  leading  edge  of  the  pulse  to  a  peak  value  at  the 
tail  end  of  the  pulse.  Behind  the  pulse,  the  density 
drops  from  the  peak  level  due  to  various  loss  mech¬ 
anisms.  The  envisioned  density  profile  is  also  plot¬ 
ted  in  fig.  1.  Therefore,  da)lc/dx<0  in  the  pulse  in¬ 
terval  and  consequently,  d(ct>2-o;^ )/d/>0.  This 
implies  that  plasma  during  the  breakdown  can  re¬ 
main  transparent  (i.e.  w >&>,*.)  to  the  pulse  as  long 
as  CDXDpe  initially.  A  frequency  autoconversion  pro¬ 
cess  leads  to  reflectionless  pulse  propagation.  More¬ 
over,  eq.  (7)  can  be  rewritten  as 

d  2  i 


an  up-shift  in  the  carrier  frequency  of  the  pulse  is 
anticipated.  In  the  following,  a  chamber  experiment 
on  wave  propagation  is  performed  to  observe  the  fre¬ 
quency  up-conversion  phenomenon  as  a  manifesta¬ 
tion  of  eq.  (7).  It  is  noted  that  it  is  desirable  to  keep 
the  distortion  of  the  pulse  due  to  plasma  generation 
small  so  that  can  be  kept  nearly  a  constant 

value  along  the  pulse’s  trajectory  thereby  achieving 
the  greatest  frequency  shift. 


3.  Experimental  results 

The  block  diagram  of  the  experimental  setup  is 
shown  in  fig.  2.  A  directional  coupler  connected  to 
the  horn  transmitting  the  input  pulse  is  used  to  mon¬ 
itor  the  spectrum  of  the  input  pulse  which  is  used  as 
the  reference.  After  passing  through  the  chamber,  the 
pulse  is  received  by  a  second  horn  placed  at  the  op¬ 
posite  side  of  the  chamber.  The  chamber  is  pumped 
to  the  lowest  pressure  at  about  0. 1 5  Torr  in  order  to 
minimize  the  electron-neutral  collision  frequency  v, 
which  is  below  5  X 1 0®  s~ 1 .  Limited  by  the  available 
power  of  the  magnetron  tube,  the  ionization  fre¬ 
quency  v j  is  less  than  1 08  s" 1 .  Thus,  v  v  «  w 2.  We  also 
limit  the  power  level  of  the  input  pulse  in  order  to 
prevent  the  overdense  plasma  being  generated  near 
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the  incident  boundary,  i.e.  cupe(at.r=0)  The  as¬ 
sumptions  of  eq.  (7)  are  than  justifiable  with  the 
current  experiment. 

One  microsecond  pulse  is  used  in  the  experiment, 
thus,  the  spatial  length  of  the  pulse  is  about  300  m 
long  which  is  much  longer  thar  the  length  0.6  m  of 
the  chamber.  Since  in  any  instance  only  a  small  frac¬ 
tion  of  the  pulse  can  fit  inside  the  chamber  where 
plasma  is  generated  by  the  field  of  the  same  pulse, 
we  can  only  anticipate  to  observe  a  very  small  per¬ 
centage  of  frequency  shift  in  the  present  experiment. 
In  other  words,  one  may  divide  the  pulse  into  500 
segment,  ind  each  segment  of  the  pulse  experiences 
a  frequency  shift  while  it  is  inside  the  chamber. 
Therefore,  the  maximum  frequency  shift  on  the  out¬ 
put  pulse  may  just  be  the  maximum  frequency  shift 
of  each  segment. 

Three  sets  of  frequency  spectra  are  presented  in 
figs.  3a-3c.  In  each  set  of  frequency  spectra,  curve  1 
is  the  frequency  spectrum  of  the  input  pulse,  while 
curve  2  is  the  frequency  spectrum  of  the  output  pulse. 
The  results  are  recorded  for  three  consecutively  in¬ 
creasing  amplitudes  of  the  input  pulse  transmitted 
into  the  chamber.  The  frequency  difference  between 
the  peaks  of  each  spectral  set  represents  the  fre¬ 
quency  shift  attributed  to  the  plasma  generation  in 
the  chamber  by  the  corresponding  input  pulse.  It  is 
shown  that  the  frequency  shift  in  the  first  two  cases 
( figs.  3a  and  3b )  is  much  larger  than  that  in  the  third 
case  (fig.  3c).  On  the  other  hand,  one  can  also  note 
that  the  distortion  of  the  pulse  shape  and  the  atten¬ 
uation  of  the  pulse  intensity  in  the  first  two  cases  are 
much  smaller  than  those  in  the  third  case.  It  suggests 
that  if  the  pulse  is  not  distorted  seriously  (i.e. 
oj>  cu^)  at  the  incident  boundary,  the  pulse  will  not 
be  distorted  seriously  and  is  experiencing  frequency 
up-shift  while  ionizing  the  air  during  the  subsequent 
propagation  period.  This  is  precisely  the  implication 
of  eq.  (7). 

It  is  noted  that  the  achievable  frequency  of  the 
present  experiment  is  much  less  than  that  reported 
in  the  previous  work  [  10  j  conducted  in  the  same  ex¬ 
perimental  setup.  The  reason  is  that  the  previous  ex¬ 
periment  is  performed  by  using  two  crossed  micro- 
wave  pulses  so  that  the  plasma  in  the  chamber  is 
generated  simultaneously  in  the  intersecting  region 
of  the  pulses.  Thus,  the  initial  effect  (causing  the 
change  of  frequency)  is  dominant  over  the  boundary 
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Fig.  3.  The  frequency  spectra  of  the  input  (curve  I )  and  output 
(curve  2)  pulses,  the  background  ait  pressure  is  0.15  Too.  (a) 
Microwave  power P~  I, IF*,  the  breakdown  threshold  power,  (b) 
F~l.27F*;(c)F~1.42FlB. 


effect  (causing  the  change  of  the  wavenumber).  On 
the  other  hand,  the  present  work  is  aimed  at  studying 
and  improving  pulse  propagation,  hence,  only  a  sin¬ 
gle  pulse  is  used  in  the  experiment.  Thus,  the  bound- 
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ary  effect  and  the  initial  effect  are  always  mixed  dur¬ 
ing  the  density  change.  The  boundary  effect  works  to 
countervalance  the  initial  effect  in  causing  the  fre¬ 
quency  shift. 


4.  Discussion 

A  pulse  propagation  experiment  for  demonstrat¬ 
ing  the  frequency  autoconversion  process  in  an  ion¬ 
izing  medium  is  conducted  in  a  large  chamber.  The 
results  show  that  the  output  pulse  with  less  distor¬ 
tion  by  the  self-generated  plasma  experiences  more 
frequency  up-shift.  It  is  consistent  with  the  impli¬ 
cation  of  eq.  (7)  which  indicates  that  if  toxo^  ini¬ 
tially  and  pulse  energy  depletion  by  the  self-gener¬ 
ated  plasma  is  small,  then  toxv^  throughout  the 
propagation.  Moreover,  to  increases  along  the  pulse’s 
trajectory. 

However,  how  are  the  experimental  results  related 
to  the  point  of  frequency  autoconversion  process  for 
reflectionless  pulse  propagation?  It  should  be  noted 
that  the  air  pressure  inside  the  chamber  is  more  or 
less  uniform.  On  the  other  hand,  the  air  pressure  from 
the  ground  to  space  decreases  with  the  altitude  and 
the  breakdown  threshold  of  the  air  is  pressure  de¬ 
pendent.  In  general,  the  breakdown  threshold  de¬ 
creases  with  altitude  before  reaching  an  altitude  hav¬ 
ing  the  Paschen  breakdown  minimum.  Beyond  that 
altitude,  the  breakdown  threshold  increases  with  the 
altitude  [11].  The.refore,  if  the  pulse  is  transmitted- 
from  the  ground  to  space  and  if  the  ionization  den¬ 
sity  on  the  ground,  where  the  breakdown  threshold 
is  very  high,  does  not  exceed  the  critical  density  set 
by  the  carrier  frequency  of  the  pulse,  the  pulse  is  ex¬ 
pected  to  propagate  through  the  Paschen  minimum 
region  without  experiencing  cutoff  reflection  even  if 
the  pulse  intensity  in  that  region  is  many  times  the 


breakdown  threshold  intensity.  This  is  because  the 
carrier  frequency  of  the  pulse  is  continuously  shifted 
upward  as  more  dense  plasma  is  generated  during 
the  propagation.  Likewise,  such  a  frequency  auto¬ 
conversion  process  also  works  for  pulse  propagation 
from  space  to  the  ground. 

The  frequency  autoconversion  process  is  expected 
to  make  the  operation  of  the  high  power  radar  and 
directed  energy  systems  easier  and  more  effective. 
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ABSTRACT 


The  absorption  and  scattering  of  microwaves  by  a 
single  column  arc  plasma  are  studied  and  compared 
with  the  theory.  The  plasma  is  produced  in  a 
background  pressure  1  to  10  Torr  through  steady  state 
hollow  cathode  discharge  where  the  cathode  is  made  of 
Tantalum  tubing  allowing  Argon  gas  to  flow  through. 
The  diameter  of  the  plasma  column  can  be  varied  from 
2  to  10  cm  by  varying  the  gas  flow.  The  current  and 
voltage  of  the  arc  are  about  10  A  and  50  V,  respectively. 
This  plasma  is  shown  to  give  rise  to  strong  absorption 
and  scattering  of  incident  microwave.  The  scattering 
experiment  is  performed  by  placing  the  vacuum 
chamber  in  the  far  field  zone  of  the  transmitting  horn 
and  measuring  the  angular  distribution  of  the  scattered 
field  intensities.  The  result  shows  the  expected 
attenuation  of  the  signal  in  the  forward  direction 
together  with  the  scattering  of  the  incident  wave  to  the 
sides. 


EXPERIMENTAL  SETUP 


p 


f 


i 


Hollow  Cathode  Discharge 

The  plasma  discharge  experimental  set-up  is 
shown  in  figure  1.  A  vacuum  system  made  of  a  large 
bell-jar  was  used  to  contain  the  argon  plasma  produced 
by  hollow  cathode  discharge.  The  hollow  electrodes  are 
made  from  Tantalum  for  its  relatively  low  work 
function  and  high  melting  point.  The  background 
pressure  can  be  varied  from  300mTorr  to  about  10  Torr 
by  regulating  the  amount  of  argon  gas  flowing  through 
the  hollow  cathode.  Figure  2  shows  the  electrical 
schematic  of  the  discharge  circuit.  When  the 
background  pressure  is  stabilized  at  about  1  Torr,  the 
application  of  an  RF  trigger  to  the  electrodes  heats  the 
cathode  sufficiently  to  allow  a  self-sustained  discharge 
to  begin.  The  electrical  operating  parameters  of  the 
plasma  were  10  amps  <  id  50  volts  established  by  the 
external  power  supply  and  circuitry  which  were 
variable  up  to  100  amps  and  300  volts  maximum. 
Using  optical  and  spectroscopic  diagnostics,  the  plasma 
parameters  were  estimated  to  be  ne  ~  10ll-K)13  cm-3 
and  Te  ~  l-2eV.  By  varying  the  gas  flow,  and  thus  the 
background  pressure,  the  plasma  can  be  varied  to 
appear  as  a  pinched,  dense  column  with  diameter  of 
the  order  of  2  cm  to  a  less  dense,  elongated  spheroid 
with  minor  axis  of  the  order  of  10  cm.  The  length  of 
the  column  in  all  cases  can  be  varied  from  5  to  10 
inches  by  adjusting  the  distance  between  the  electrodes. 
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Figure  1.  Schematic  diagram  of  single  hollow  cathode  discharge  unit 
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rf  trigger :  short  to  dc 


Fig.  2  Schematic  diagram  of  the  HCD  discharge  circuit. 


Microwave  Scattering 


The  microwave  scattering  set-up  is  shown  in  figure 
3.  An  HP  microwave  sweep  generator  was  set  up  for  10 
dBm  cw  nominal  power  output  at  frequencies 
selectable  from  3.8  to  12  GHz.  The  signal  is  fed  to  a 
launching  horn  aimed  at  the  plasma  chamber  and  is 
placed  near  it.  Microwave  absorbers  are  then  placed 
surrounding  the.  chamber  to  minimize  unwanted 
reflections  and  scattering  from  nearby  objects. 
Windows  in  the  absorbers  allow  the  placement  of  a 
receiving  horn  at  various  angles  about  the  center  of  the 
plasma  chamber  to  pick  up  the  microwave  signal.  A 
spectrum  analyzer  also  capable  of  time  domain  sweeps 
is  then  used  to  examine  the  received  signal  in  the 
presence  and  in  the  absence  of  the  plasma  column. 
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Fig.  3  Schematic  diagram  of  experimental  setup. 


EXPERIMENTAL  RESULTS 


Broadband  Attenuation 

A  typical  frequency  domain  sweep  from  4-8  GHz 
comparing  the  forward  received  signal  with  and 
without  plasma  is  shown  in  figure  4.  It  is  seen  that 
attenuation  to  the  transmitted  signal  in  the  forward 
direction  attributable  to  absorption  and  scattering 
occurs  over  the  whole  frequency  range.  Fluctuations  in 
the  amount  of  attenuation  are  due  to  fluctuations  in 
the  plasma  during  the  time  the  spectrum  analyzer  took 
its  reading  at  each  particular  frequency.  This  plasma 
fluctuation  has  been  observed  experimentally  by  the 
use  of  a  photomultiplier  tube  and  the  cause  has  been 
traced  to  the  ac  ripples  in  the  source  supply.  Since  the 
sampling  algorithm  in  the  analyzer  displays  the  peak 
received  value  within  any  sampling  period,  one  expects 
that  the  maximum  observed  attenuation  should  apply 
to  all  frequencies  if  the  plasma  can  be  made  more 
quiescent. 
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Figure  4  Spectral  distribution  of  transmitted  wave.  (4-8  GHz  sweep) 
Curve  1  -  without  plasma  column  in  chamber. 

Curve  2  -  plasma  column  present. 


Angular  Distribution  of  Scattering  Wave 


Measurements  of  the  angular  spectrum  of 
microwave  scattering  by  the  plasma  column  at  6  GHz 
are  shown  in  figure  5.  Three  sets  of  data  taken  under  as 
close  as  possible  conditions  were  used  to  generate  the 
curve  and  error  bars.  The  data  points  are  normalized  to 
the  incident  wave  for  comparison  with  the  theoretical 
results  to  be  presented  in  the  future.  It  can  be  seen 
from  the  attenuation  of  the  signal  in  the  forward 
direction  and  the  enhancement  of  received  signal 
above  20  degrees  from  forward  that  the  plasma  column 
indeed  scatters  the  incident  wave  away  from  the 
forward  and  to  the  sides.  The  attenuation  to  the 
incident  wave  in  the  forward  direction  as  shown  is  26 
dBs.  The  typical  value  of  the  enhanced  signal  is  -10  dB 
below  the  incident  power  level. 
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Figure  5  Angular  distribution  of  scattering  wave 


